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Abstract A number of different families of nanowires 
which self-assemble on semiconductor surfaces have been 
identified in recent years. They are particularly interest- 
ing from the standpoint of nanoelectronics, which seeks 
non-lithographic ways of creating interconnects at the 
nanometre scale (though possibly for carrying signal rather 
than current), as well as from the standpoint of tradi- 
tional materials science and surface science. We survey 
these families and consider their physical and electronic 
structure, as well as their formation and reactivity. Par- 
ticular attention is paid to rare earth nanowires and the 
Bi nanoline, both of which self-assemble on Si(OOl). 



1 Introduction 

As the scale of architectures in integrated circuit design 
continues to be reduced, the dimensions of present-day 
interconnects are on the scale of tens of nanometers, and 
the interconnects for next-generation nanoelectronic de- 
vices may well be just a few nanometers in lateral dimen- 
sion. Furthermore, the incorporation of nanometer-scale 
components, such as active molecules, into integrated 
circuits will require interconnects of a similar scale. On 
these scales, self-assembled nanowire systems are becom- 
ing increasingly interesting^|2], as conventional litho- 
graphic techniques reach their limits around 10-15 nm 3 
and SPM-based nanolithography|3j methods lack scala- 
bility. Moreover, the introduction of such "bottom-up" 
technology, based on naturally nanometer-scale compo- 
nents to complement or even replace the current "top- 
down" technology is likely to require completely different 
architectures. One example is the "crossbar" architecture jH], 
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in which active molecules are used as devices at the junc- 
tions between two perpendicular nanowires. This archi- 
tecture is designed to take advantage of the typical prod- 
uct of self-assembly schemes - an array of parallel wires - 
rather than relying on controlled positioning of individ- 
ual wires. Over the last ten to fifteen years, a number of 
systems have been identified which yield self-assembled 
nanowires on semiconductor surfaces, and stem from a 
mixture of serendipity and deliberate design. There are 
only a few systems which provide long, straight nano- 
lines on flat terraces. In/Si(lll)|n] (amongst other met- 
als) and Pt/Ge(001)|^ are examples of systems which 
undergo surface reconstructions leading to ID atomic 
chains with metallic properties. On the Si(OOl) surface, 
the Bi nanolines|H||n] appear to be the only example of 
an atomically-perfect, straight, self-assembled ID nanos- 
tructure, but are not metallic. Aside from fortuitous dis- 
coveries, symmetry-breaking by use of a vicinal surface 
to produce and array of steps can be used to grow long, 
straight nanowires of many metals, including Au, Bi et 
al.^Ol, and 1-D nanoscale structures can be designed 
by taking advantage of anisotropy in the heteroepitax- 
ial strain between the Si(OOl) surface and an appropri- 
ate deposited material. The rare-earth silicidcs arc the 
archetype of this method^J. Within the family of rare- 
earths (as well as scandium and yttrium) , a close lattice 
match can be obtained in one direction, with a variety of 
positive and negative lattice mismatchs in the orthogo- 
nal direction. By choosing the right material, metastable 
ID nanoscale structures 3-10 nm wide and hundreds of 
nanometres long can be formed. These have some phe- 
nomenological similarities with the Bi nanolines, as will 
be described. This review will survey the growing field 
of nanowires on semiconductor surfaces, with the focus 
on the Bi/Si(001) nanoline system, whose importance is 
explained in the next section. 

The study of semiconductor surfaces has been revolu- 
tionised over the last twenty years by the advent of two 
techniques which give real-space atomic scale informa- 
tion: Scanning Tunneling Microscopy (STM) and elec- 
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tronic structure modelling (in particular, Density Func- 
tional Theory, or DFT) . The two techniques are comple- 
mentary, in that STM gives an approximate answer to 
an exact question (i.e. the structure of the real Si(OOl) 
surface), while DFT gives an exact answer to an ap- 
proximate problem (i.e. the lowest-energy structure of a 
computation cell containing a few Si dimers, and a few 
layers of bulk-like Si). While these techniques between 
them enabled the clear identification of the Si(OOl) re- 
construction, even today there is some controversy about 
the exact nature of the buckling of dimers on Si(OOl) 
at low temperatures IT2I. In the same way, the under- 
standing of the structure and properties of the nanolines 
surveyed in this paper has been greatly enhanced by a 
close collaboration between experiment and theory. The 
benefits of this synergistic approach will be drawn out 
throughout the discussion. 

1.1 Introducing the Bi nanolines 

The Bi:Si(001) nanoline system, which will be the focus 
of this review, was discovered about ten years ago. As 
can be seen from the examples in Fig. ^ they grow per- 
fectly straight along the (110) directions on the Si(OOl) 
surface for hundreds of nanometres, apparently limited 
only by the terrace size of the underlying substrate. 
When they encounter a step edge, they will either grow 
out over lower terraces in a long, narrow promontory, or 
they will burrow into higher terraces, until a deep inlet 
is formed. They have a constant width, 1.5 nm or four 
substrate dimers, and are very stable. So long as the tem- 
perature is not so high that the Bi can evaporate, they 
are stable against prolonged annealing, maintaining the 
same width. This is quite unlike the rare earth silicide 
systems which in most cases (Yttrium being the proba- 
ble exception) will coarsen into 3D epitaxial islands with 
long anneals. There are no other elements which produce 
perfect nanolines on the Si(OOl) surface, although there 
is some evidence to suggest that Sb nanowires with the 
same appearance as the Bi nanoline form in the pres- 
ence of surface hydrogen jl^, and DFT modelling that 
indicates that a Sb nanoline would be stable J3j. There 
has also been a suggestion that Er may form nanolines 
of the same structure 15 (as opposed to silicide wires), 
but this remains unconfirmed. 

Structurally, the Bi nanoline is unique amongst nano- 
line systems — it is neither a periodic reconstruction of 
the surface, nor the result of anisotropically strained het- 
eroepitaxial growth of one bulk structure on another. In 
fact, it is somewhere in between. The nanoline structure 
is built around a pair of Bi ad-dimers, on top of a com- 
plex subsurface reconstruction which is responsible for 
many of the nanoline's remarkable properties, such as 
the extreme straightness; this will be explained in detail 
below. In this article, we intend to describe the forma- 
tion, structure and properties of this system, setting it 




Fig. 1 A typical high-temperature STM image of the 
Bi/Si(001) surface. Here Ifimxlfim. The Bi nanolines run 
along (110) directions, perpendicular to the Si dimer rows, 
and can grow to over l^m in length. In the inset, the atomic 
structure, comprising a pair of Bi dimers, about 6.5A apart, 
is shown. Hydrogen termination of the Bi nanolines has been 
used to enhance the apparent height of the nanoline. 

in the context of other nanoline systems, particularly the 
rare-earth silicide family. 

1.2 Overview 

The paper is laid out as follows: in Section |21 we con- 
sider the methods used to examine the systems, both 
experimental and theoretical; in Section|31self-assembled 
nanowire systems other than the Bi:Si(001) nanowire 
system are described, with particular reference to the 
rare-earth silicide family; in Section ^ the Bi:Si(001) 
nanowire system is examined in detail, and in Section |S1 
the state of the field is considered and conclusions are 
drawn. In considering the nanowires in Sections |31 and 0] 
we consider first the formation and structure of the sys- 
tems (including the effects of annealing) , followed by the 
electronic properties (including conductance data where 
available) and finish with the reactivity of the systems. 

2 Methods 

2.1 Hot STM 

Since its invention in 1982 by Binnig and Rohrer at 
IBM, STM has allowed researchers to image dynamic 
surface processes, in real-space, with near-atomic resolu- 
tion. STM as a technique bridges a variety of disciplines. 
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In surface physics, STM is unique as a real-space probe 
of electronic density of states of individual species; the 
STM tip can be a passive observer of surface chemical 
reactions, and also an active catalyst of these reactions; 
and STM allows for the direct observation of many phase 
transformation processes, albeit in 2D rather than 3D. 
However, this is not meant to be a review of capabil- 
ities of STM; for that, the reader is pointed to classic 
books on the subiect [TCIll7| . Here a few comments are 
made, which are pertinent to the particular challenges of 
elevated-temperature ('hot') STM for the study of epi- 
taxial growth. 

Surface reactions can be followed at room temper- 
ature by using quench experiments, in which a high- 
temperature surface is cooled rapidly to freeze in a snap- 
shot of an evolving surface, for example in the growth of 
HI-V semiconductor surfaces ^Bl- However, the quench- 
ing process may introduce spurious surface features, and 
so it is preferable to image the surface in-situ, i.e. at tem- 
perature. The ability to vary the sample temperature 
makes it possible to image a surface in the process of 
changing, rather than after it has changed. To give just 
one example, the phase transition between the (1x1) re- 
construction and the (7x7) reconstruction on the Si(lll) 
surface, at around HOOK, has been observed directly 
in both directions, by careful control of the tempera- 
ture across the transition 19^ For useful observations of 
dynamic processes, the speed of the reaction must be 
matched to the rather slow imaging rate of an STM. 
One example is the motion of step edges [30]. At around 
300°C, step motion can be imaged as the addition or 
subtraction of 4 atoms between STM images. However, 
above about 500°C, the step edge will have moved sig- 
nificantly in the time taken for the STM to scan one line. 
Thus from line to line, the step will be in different places, 
as it oscillates about a mean position. In STM, there- 
fore, fast-moving steps have a streaky appearance above 
500°C, as may be seen in Fig. El Furthermore, mobile 
species may be invisible at elevated temperatures, either 
due to the decrease in contrast, or due to their mobility. 

Surface chemical processes, such as the diffusion of Si 
atoms on Si(OOl) 20 or the decomposition of a molecule|21j. 
can be studied at the single-species level, and complete 
chemical reaction pathways from initial adsorption can 
be followed through to complete decomposition |221ES|- 
Kinetic information such as activation barriers can be 
extracted, and any metastable intermediate structures 
can be imaged, giving a unique insight into the reac- 
tion process. The ability to control the substrate tem- 
perature is particularly important in studies of epitaxial 
growth, where the kinetic pathways available at differ- 
ent temperatures will determine the surface morphology. 
In the case of Si/Si(001) and Ge/Si(001) epitaxy, it has 
been possible to image the surface simultaneously with 
the deposition of material, and therefore observe directly 
the nucleation of islands, and the growth of these nuclei 
into ID, 2D and 3D islands, along with associated pro- 



cesses, such as the ripening of epitaxial islands, surface 
stress relief mechanisms, and the transition from island 
growth to step flow growth modes |i24«25u26) . all of which 
processes occur at different sample temperatures. 

However, these added abilities bring with them a 
penalty in sensitivity. There is generally a loss of res- 
olution due to thermal noise above room temperature, 
and a loss of stability, with thermal drift making it hard 
to image the same area for extended periods. In many of 
the high-temperature STM images shown here, the nano- 
lines have a curved appearance. They have no curvature 
themselves, the appearance comes from the drift dur- 
ing scanning. Furthermore, imaging of epitaxial growth 
necessarily involves imaging a surface where there are 
many atoms loosely bound to the surface, or in the gas 
phase being deposited. In this situation, it is very easy 
for atoms to stick to the STM tip, which can have ad- 
verse effects on the imaging, and it is usually prefer- 
able to obtain stable imaging conditions rather than 
achieve the highest possible resolution, which is generally 
a more unstable situation. For the same reason. Scan- 
ning Tunneling Spectroscopy (STS), or Current Imag- 
ing Tunneling Spectroscopy(CITS)^, becomes unfeasi- 
ble at elevated temperatures. It is generally better to 
take a series of images of the same area at different 
bias voltages, giving voltage contrast of an object, al- 
though as the feedback loop is active in this situation, 
the tip-sample distance will vary, and so the information 
gained is not equivalent to STS/CITS. Thus these limita- 
tions mean that while high-temperature STM brings the 
unique ability to study evolving surface morphology in 
real-space, it is not always the best way to study growth 
surfaces. A combination of high-temperature studies, and 
room-temperature studies of quenched surfaces, is likely 
to provide a fuller description of a system. 

2.2 Electronic Structure Modelling 

There have been a number of books on modelling tech- 
niques published recentlv j?7ll28| . so this section will fea- 
ture only a brief description of the techniques used for 
modelling of semiconductor surfaces, concentrating on 
electronic structure calculations (both semi-empirical and 
ab initio) . These techniques retain the quantum mechan- 
ics of the electrons while treating the ions classically. 

The first problem to consider is how to model a semi- 
infinite piece of material (that is, a surface). The approx- 
imation that is used is to consider a small piece of ap- 
propriate material with boundary conditions. Generally, 
one of two solutions is used: either a cluster of atoms, 
with the dangling bonds on the edges terminated in some 
suitable way (e.g. hydrogen), or a slab of atoms with pe- 
riodic boundary conditions in two or three directions. A 

^ A process in which a rapid voltage ramp is applied to the 
STM tip while in a fixed position above the surface, gener- 
ating an I/V curve from which the LDOS can be measured 
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cluster allows a smaller number of atoms to be used than 
a slab, but has the drawback that it is hard to allow for 
long-range strain effects (such as those caused by recon- 
structions or steps on surfaces). A slab approach requires 
care: the vacuum gap between repeating images must be 
large enough to prevent interactions, and the slab itself 
must be sufficiently thick to allow strain relaxation. The 
slab approach is more commonly used, particularly with 
the methods described below. 

Once the system to be modelled has been defined, 
the calculation technique must be chosen. Putting aside 
quantum chemical methods for the sake of brevity (they 
scale rather poorly with system size — for more details, 
the reader is directed to one of many excellent books 
on the subject[2ni), we will consider the semi-empirical 
tight binding (TB) technique, and the ab initio density 
functional theory (DFT) technique. 

TB .UTITTT] postulates a basis set of atom-centred or- 
bitals^ for the wavefunctions; the Schrodinger equation 
can then be rewritten as a matrix equation, with ele- 
ments of the Hamiltonian matrix formed from integrals 
between orbitals on different atoms. These Hamiltonian 
matrix elements are fitted to either ab initio calculated 
data or experimental results. For simplicity, the basis 
is often assumed to be nearest neighbour only (with a 
cutoff defined on the range of interactions) and orthog- 
onal. Once the Hamiltonian has been defined, the band 
energy of the system can be found by diagonalisation 
of the Hamiltonian or other methods. The other ener- 
getic terms (e.g. Hartree correction, exchange and ion- 
ion interactions) are represented by a repulsive potential 
(often a pair potential). 

Despite its apparent simplicity, TB is extremely ef- 
fective, and generally qualitatively accurate, if not ap- 
proaching quantitative. It has been shown[31] that this 
TB formalism can be derived from DFT via a set of 
well-defined approximations, which explains to some ex- 
tent its success. It is well-suited for a rapid exploration 
of configuration space (e.g. possible structures for some 
new surface feature jSS]) provided that a suitable param- 
eterisation exists for the bonds between different species. 
The fitting of parameterisations is non-trivial, and the 
transferrability of a given parameterisation (i.e. its accu- 
racy in environments far from those in which the fitting 
was performed) is never guaranteed. Generally the ma- 
trix elements themselves are fitted to a band structure 
(or energy levels of a molecule), while their scaling with 
distance is fitted to elastic constants or normal modes 
of the system. TB retains quantum mechanics (since the 
energy is obtained by solving the Schodinger equation) 
while approximating the most complex problems. 

DFT[nSlEZ| starts with an exact reformulation of the 
quantum mechanics of a system of interacting electrons 
in an external potential: the result is a set of equations 

^ For the class of methods known as ab initio tight 
binding [32i33f . the basis is explicit. 



for non-interacting electrons moving in an effective po- 
tential with all the complex electron-electron interac- 
tions in a single term, known as the exchange-correlation 
functional (it is a function of the charge density, which 
is itself a function of position, hence "functional"). Un- 
fortunately, the form of this functional is not known, 
and must be approximated, for instance with the local 
density approximation (LDA) or one of the generalised 
gradient approximations (GGA). 

DFT has been extremely successful in many areas of 
physics, chemistry, materials and, increasingly, biochem- 
istry, particularly when combined with the pseudopoten- 
tial approximation. The hard nuclear potential and the 
core electrons of each atom are replaced with a single 
"pseudopotential" , and only the valence electrons are 
considered, leading to a softer potential. This approxi- 
mation is most effective for atoms where there is con- 
siderable screening of the valence electrons by the core 
electrons (e.g. in Si the 3s electrons are screened by the Is 
and 2s shells, while the 3p electrons are screened by the 
2p shell); in first row elements and first row transition 
metals, this is a much smaller effect, leading to the de- 
velopment of "ultrasoft" pseudopotentials. These issues 
are discussed in more detail in many excellent books and 
reviews P51I77| . The essential point to note is that DFT 
is widely used for electronic structure calculations, and 
provides results which are accurate to within approxi- 
mately 0.1 eV. 

An important development in electronic structure 
techniques over the last ten years is that of linear scaling 
techniauesP^. Standard techniques for solving for the 
ground state, whether TB or DFT, scale with the cube 
of the number of atoms in the system (either because 
matrix diagonalisation scales with the cube of the ma- 
trix size, or because the eigenfunctions spread over the 
whole simulation cell, leading to cubic scaling when they 
are orthonormalised) . This scaling places rather strong 
restrictions on the sizes of system which can be mod- 
elled, even on massively parallel machines: for ab ini- 
tio methods, going beyond 1,000 atoms rapidly becomes 
prohibitive, though somewhat larger systems can be ad- 
dressed. However, since electronic structure is fundamen- 
tally local (consider bonding as an example), the amount 
of information in the system should be proportional to 
the number of atoms in the system. Tight binding meth- 
ods which take advantage of this (e.g. refs. j4()ll41| 'l have 
been widely used for some time, allowing calculations on 
many thousands of atoms. The implementation of linear 
scaling DFT algorithms has proved significantly harder, 
though recent efforts in nearly linearly and linear scal- 
ing methods 43 4 jll45ll4f)| suggest that these techniques 
are starting to produce useful and general results. 

2.3 Modelling STM 

One of the challenges of using STM to examine semicon- 
ductor surfaces is that the current arises from both geo- 
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metric and electronic structure, though this is less true 
at high biases (where height or geometric structure will 
dominate). Some technique is required to understand 
these changes, and to test proposed structures against 
experiment. The field of STM simulation is a complex 
one (in part because the structure and composition of 
the tip is unknown); the interested reader is referred 
to excellent reviews for further information j47ll48| . We 
will briefly summarise the simplest and most common 
approximation used, and discuss how the experimental- 
theoretical interaction can best proceed. 

The most commonly used approach is the Tersoff- 
Hamann approximation^^!- This asserts (via a careful 
series of approximations) that the tunneling current is 
proportional to the local density of states (LDOS) due 
to the sample at the position of the tip. In effect, the tip 
is assumed to have a flat density of states. This approach 
is directly equivalent to considering the projected charge 
density (i.e. the partial charge density due to each band 
for bands within Vbias of the Fermi level). It is qualita- 
tively accurate, though will not reproduce the observed 
corrugation values correctly |48j. 

The authors are of the opinion that for true success 
in investigating nanowires (and other systems) on semi- 
conductor surfaces, a close collaboration between exper- 
iment and modelling is required. Any successful collab- 
oration depends on a number of factors, but the key 
factors, we believe, are frequent, easy correspondence, 
an understanding of the limitations and abilities of the 
appropriate techniques, and trust between the differ- 
ent sides in the collaboration. The last point extends to 
sharing of unpublished data, and extensive discussion of 
possible courses of action. The three authors have been 
working together in different combinations for ten years 
(leading to around 20 joint publications), and we find 
that our collaboration is still immensely fruitful. 

3 Other Nanowires 

While the Bi/Si(001) nanolines have remarkable struc- 
tural qualities, they are not the only self- assembled nanowire 
system on semiconductor surfaces. In this section we 
present an overview of these other nanowire systems 
(more details are presented elsewherejj), and we will 
consider in more detail the most important of these sys- 
tems, the rare-earth silicidc family, focussing on their 
structure and formation, electronic properties and reac- 
tivity. 

A variety of metals have been found to self-assemble 
into nanowires or nanolines on the Si(OOl) surface, though 
none with the perfection of the Bi nanolines. The sim- 
plest way in which to form ID structures on the Si(OOl) 
surface is by epitaxial growth. An island nucleus is essen- 
tially an ad-dimer, which sits between the dimer rows. 
The two pairs of Si dimers which support the ad-dimer 
are distorted by its presence, and so are attractive ad- 
sorption sites for further ad- atoms or ad-dimers|Sn|- This 




Fig. 2 (a): A dangling-bond wire on the H:Si(001) surface. 

(b) Ga deposited along a row of dangling bonds on H:Si(001). 

(c) : A line of Si dimers after removal of H termination by an 
STM tip. (d): Preferential adsorption of Ag onto the clean 
Si dimers, to form an atomic-scale wire. (a,b) courtesy of T. 
Hashizume, Hitachi ARL, Hatoyama, Saitama, Japan. (c,d) 
courtesy of T.Sakurai, NIMS, Tsukuba, Ibaraki, Japan. 



process has been described as a surface polymerisation 
reactionist. The long sides of the string, which are equiv- 
alent to A-type steps on this surface, have a very low 
sticking coefficient which keeps the island from broaden- 
ing at lower growth temperatures. Deposition of Group 
III and Group IV elements on Si(OOl) will therefore re- 
sult in long ID chains of dimers |S2ll53ll54ll55ll56j . How- 
ever, there is little control over the length or structure 
of these wires, and they are unstable against anneal- 
ing, eventually reorganising into compact islands. The 
Pt/Ge(001) system|7j forms nanowires through a com- 
plex surface reconstruction, resulting in a system with 
a high degree of perfection, which approaches that of 
the Bi nanolines. When annealed to high temperatures, 
arrays of atomically perfect nanowires of Pt/Ge form, 
each of which is 0.4 nm wide with a spacing of 1.6 nm 
between the wires. 

Sideways growth of two-dimensional islands can also 
be blocked by introduction of missing dimer trenches. 
These form to relieve surface stress, cither because of 
contamination by a small amount of a transition metal 
such as Ni, or through heteroepitaxial growth of another 
material such as Ge. The trenches tend to line up in 
a semiregular array, giving an approximate (2xn) re- 
construction in LEED, with n'^8-12 (though the regu- 
larity of the trenches is rather poor). This method has 
been used to deposit a variety of metals, such as Fe|57|. 
Ga|SHI and Inj59j (using the Ni-based technique) and 
the molecule styrenelHS]. The result is long-range, rea- 
sonably well ordered wires of the metals, though their 
properties have not been characterised. 
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The 2D symmetry of the surface can also be broken 
deliberately by the use of vicinal substrates to produce 
regular arrays of stepsjfil) which can then be decorated 
with materials such as gold |lUII62ll63ll64ll65| and Bi|^. 
This approach is used primarily on Si(lll), although 
with an extreme orientation towards (112) directions the 
surface might be labelled Si(557) or Si(5 5 12) depending 
on the angle. The resulting step reconstructions comprise 
chains which are metallic and one-dimensional; early ob- 
servations suggested possible observations of Luttinger 
liquid behaviour 1671, though this has been disputed |62l 
Ifi4j . Their spacing can be controlled to some extent by 
varying the angle of miscut of the surface. 

As an atomic-scale extension of the idea of lithogra- 
phy, the Si(OOl) surface has been passivated with atomic 
hydrogen (forming the monohydride phase) and individ- 
ual hydrogen atoms removed with an STM tip to form 
atomic-scale patterns. The wire formed by removal of 
hydrogen atoms along (or across) a dimer row is known 
as a "dangling bond" wire (or DB wire)'fe8 %OII70ll711 
I72II73| and is predicted to show conduction effects simi- 
lar to conjugated polymers with polaronic and solitonic 
effects |7^75ll7fij . An example of a dangling-bond wire is 
shown in Fig. Efa) . These wires have been reacted with 
a variety of adsorbates including iron [77], gallium |lir75| . 
aluminium'79^ silver "SI? and organic molecules such as 
norbornadiene 81,82,83, . However, although for a sin- 
gle row of dangling bonds, a perfect atomic wire can be 
formed, as with Ga in Fig. Efb), for greater widths, as 
in Fig. |2c,d), the dangling bond wire is more ragged, 
and the resulting wires (here Ag) are typically com- 
posed of a large density of small, roughly spherical crys- 
tals, which exhibit many imperfections and boundaries. 
A more elegant implementation stems from careful ap- 
plication of molecular chemistry to produce a directed 
reaction: when the monohydride surface was exposed to 
styrene and one hydrogen atom removed, self-directed 
growth of lines of styrene resulted from a chain reac- 
tion between each adsorbed molecule and an adjacent 
hydrogen|H2- Such self-organised methods are important 
for the growth of molecular nanowires on surfaces. How- 
ever, for all these systems there is a lack of control over 
the size and length of the nanowires. 

3. 1 Structure and Formation of Rare-Earth Nanowires 

A more promising materials system for nanowire growth 
on Si(OOl) is that of the rare-earth metal silicides. There 
has been extensive research on these systems because of 
their good conductivity and low Schottky barrier with 
siliconjSS], though this work concentrated on the Si(lll) 
surface. Moreover, in the hexagonal phase some rare- 
earth silicides show a good lattice match with the Si(OOl) 
substrate in one direction, with a large mismatch in the 
other. Heteroepitaxial growth of these materials would 
therefore be expected to be constrained in the high mis- 
match direction, and facile in the other, resulting in long 
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Fig. 3 A large-scale image of Dy NWs, showing the 3 differ- 
ent types of surface feature associated with the RENWs. On 
the terraces, a reconstruction is typically seen, here (2x7). 
Single NWs grow in both (110) directions, blocking each 
other's growth. Finally, the NWs can form bunches. Both the 
single NWs and the bunches can grow a second layer of mate- 
rial, unlike the Bi nanolines. Image courtesy Prof. J.Nogami, 
U. Toronto, Canada. 



Rare earth 


a(A) (%) 


c (A) (%) 


ScSii.7 

YSi2 (hexagonal) 
SmsSis (hexagonal) 
GdSi2 (hexagonal) 
DySi2 (hexagonal) 
HoSi2 (hexagonal) 
ErSi2-a: (hexagonal) 
YbSi2 


3.66 (-4.69) 
3.842 (+0.05) 
3.90 (+1.64) 
3.877 (+0.96) 
3.831 (-0.23) 
3.816 (-0.63) 
3.79 (-1.30) 
3.784 (-1.46) 


3.87 (+0.78) 
4.144 (+7.92) 
4.21 (+9.64) 
4.172 (+8.65) 
4.121 (+7.32) 
4.107 (+6.95) 
4.09 (+6.51) 
4.098 (+6.71) 



Table 1 Lattice mismatch between hexagonal phase of var- 
ious rare-earth silicides and Si(OOl). Data is presented in 
terms of lattice constant and percentage mismatch (the 
Si(OOl) surface has a lattice constant of 3.84A). Where the 
indication "hexagonal" is given, it indicates that other phases 
are possible. Courtesy of Dr C. Otibuchi, NIMS, .Japan.. 



ID islands. This is indeed what happens. The behaviour 
of the rare-earth nanowire systems on Si(OOl) has been 
investigated extensively, and nanowires have been ob- 
served under certain growth conditions. As the other 
major family of self-assembled nanowires on the Si(OOl) 
surface, we give a fuller description of this family for 
comparison to the Bi nanolines. 

The most common form of rare-earth nanowire (RENW) 
formed on Si(OOl) is thought to result from the anisotropic 
strain between the AIB2 crystal structure of the nanowire 
and the Si(OOl) substrate, leading to fast growth along 
the less-strained direction, and extremely limited growth 
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along the more-strained direction. The mismatch for var- 
ious silicides in the AIB2 hexagonal structure is given 
in Table ^ RENWs of this type have been reported 

FT, 98 , 99 , 100 , lO^^'^^^^ET^^^I^a^'^'^ 
[103 . For Sm 9 0..92) NWs only appear on vicinal sub- 
strates; other reDorts [Mlll04| find that on flat surfaces 
only rectangular, 3D islands are formed after a (2x3) 
layer). Bundles of NWs have been reported for Yb jl(J5| 
; again, early reDorts [Mlll06| suggested that NWs were 
not formed, and 3D islands resulted from annealing. Thoug] 
they are not rare-earth metals, similar NWs have been 
reported for both Sc ,^9 and Y |MlllU7| . resulting from 
the same combination of crystal structure and anisotropic 
strain. Although it forms a similar wetting layer on the 
substrate to Yb, Eu does not form NWs |108j . 

The resulting RENWs are found both individually 
and in bundles; the individual wires have sizes which 
vary from element to element but are in general 5-10 nm 
wide and less than 1 nm high. The wires grow extremely 
fast, reaching lengths of up to 1/im; more interestingly 
they always run along (110) directions and will grow out 
over step edges, drawing terraces with them, features 
that they share with Bi nanolines. The large-scale im- 
age in Fig. 13 is very similar to that of Bi nanolines at 
the same scale. For almost all materials, the surface of 
the wires shows a c(2x2) reconstruction 87,88. 10 3il97l 
I99lll01j . though the Yb NWs, which are off-axis relative 
to the other RENWs, may have a (1x1) surface [TT)??] . 

The Si(lll) surface matches the lattice of the RE 
AIB2 structure reasonably well in all directions (with 
a mismatch of less than 2% for most metals) 85 and 
so would not be expected to give anisotropic growth, 
but GdSi2 NWs have been formed along step edges of a 
vicinal surface [llOI . These form because of a mismatch 
perpendicular to the step edge, and grow to over 1/im 
long on appropriately prepared samples, with a width of 
10 nm and height of 0.6 nm. 

Other NW structures have been reported for tran- 
sition metals on Si(lll), for instance Fe |lllj (though 
these are rather short wires), NifTn) and TifTTSIIIl- 
The growth mechanism for these wires is not entirely 
clear: the resulting NWs are 5 nm high and 5 nm wide 
(NiSi2), 10 nm high and 40 nm wide (TiSi2), though for 
both of these systems the NWs appear to grow into the 
substrate somewhat. There is a technique for growing 
NWs on Si(OOl), Si(llO) and Si(lll) explicitly relying 
on growth into the substrate, resulting in "endotaxial" 
NWs, which has been applied to DyjnH] in the hexago- 
nal AIB2 structure, as well as Co lll5lll00| . which takes 
the CaF2 structure. The driving force appears to be ki- 
netic in these systems: the long direction of the NW has 
an interface which grows faster than the short direction, 
though this is partly dependent on the structure of the 
interface below the surface. 

The growth and formation of the RENWs is not un- 
derstood at an atomic level: they form extremely quickly 




Fig. 4 While a wetting layer is common for Dy silicide wires, 
it is not always present, and depends upon the total cover- 
age, and annealing condition. In (a), there is a reconstructed 
wetting layer, and most NW are isolated and disconnected, 
with a consistent width, although second-layer growth has 
occurred along some NWs. In (b), the background is disor- 
dered Si(2xl), and most wires are bundles, which have joined 
together where they meet, forming multilayer islands. Images 
courtesy Prof. J.Nogami, U. Toronto, Canada. 
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(on experimental timescales) so that, for instance, hot 
STM cannot be used to observe formation of the hnes. 
UnHke the Bi nanohnes, which have a constant width, 
and do not change in width with anneahng, the rare- 
earth nanohne family can take a variety of widths de- 
pending on the growth conditions. In the case of DySi2, 
single NWs of a consistent width form from the under- 
lying (2x7) reconstructed wetting layer, as in Fig. 01 
However, further annealing results in the formation of 
bunches of nanowires, which join together where they 
meet, and multilayer islands form. As strained coherent 
heteroepitaxial islands, the RENWs are only metastable; 
prolonged annealing will cause them to coarsen, and 
form large 3D islands. Overall, the deposition/anneal 
temperature must be in the window between 550°C and 
650°C, and annealing time should not exceed 2-10 min- 
utes. Below 550°C or 2 minutes, the reaction will be in- 
complete and the NWs formed will be immature. Above 
650° C or 10 minutes the systems tend towards lower 
energy, thermodynamically stable states. In particular, 
Dy-Si(OOl) will form 3D islands if annealed for long times 
(e.g. 30 minutes for 0.86ML)|nn| with wires formed for 
shorter annealing periods. Annealing Er-Si(OOl) initially 
results in NWs, but as the annealing time is increased 
dislocations form in the NWs|89|. If the temperature is 
raised beyond 620° C the same behaviour is seen, lead- 
ing ultimately to coarsening into islands E3- ^ large- 
scale PEEM study of ErSi2[nil found that large NWs 
(which may have been bundles of nanowires — the tech- 
nique does not have the resolution to distinguish) were 
not affected by annealing on a coarse scale. However, 
Gd shows more complex behaviour. Annealing studies 
of thin films of Gd found that both hexagonal and or- 
thorhombic phases could form^^, and that with longer 
annealing times the orthorhombic phase grew at the 
expense of the hexagonal phase. Long anneals of the 
Gd/Si(001) system, up to 620°C for one hour, resulted 
in the formation of previously unknown silicide struc- 
tures, which are aligned perpendicular to hexagonal sili- 
cide NWs on the same terrace [TT7) . 

There are various useful observations of their be- 
haviour which can be made: 

— A wetting layer with a characteristic reconstruction 
forms before the nanowires, and persists on the sub- 
strate with the nanowires; specific reconstructions 
have been observed for: 

- Eu-Si(001):(2x3)[Tni, 

- Sm-Si(001):(2x3)[TniI, 

- Yb-Si(001):(2x3)inH||IIl||IIllIinHl and (2x4)[TIl| 

- Nd-Si(001):(2x3) and (2x4)[Tni, 

- Dy-Si(001):(2x4)[ni, and (2x5) and (2x7) 97 , 

- Gd-Si(001):(2x4)[S3, and (2x5) and (2x7)^ 

- Ho-Si(001):(2x4) and (2x7)[TIS! 



— The only system for which this does not happen is 
Er-Si(OOl) 87 . 

— The system must be heated during deposition or an- 
nealed post-deposition to allow reaction of the rare- 
earth metal with the substrate ■86",'94 

— Annealing the system for too long, or depositing at 
too high a temperature, can result in 3D islands rather 
than nanowires; this is discussed in more detail be- 
low. 

As different metals have different properties for each 
of these categories, we will discuss them briefly in turn. 
Gd and Dy[57j form a (2x4) reconstruction first, fol- 
lowed by a (2x7) reconstruction |97,Ji}2, just prior to 
nanoline formation. Rather similar behaviour is seem for 
Ho jlOB) . The metal fraction in the two reconstructions 
is 3/8 and 5/14 ML respectively [111113 ■ For these three 
metals, the presence of the wetting layer is intimately 
connected with formation and stability of the nanowires. 

Despite the wealth of studies of ErSi2 nanowires, 
there are no reports of the structure and make-up of a 
wetting layer for this system; there is some evidence j87) 
that once more than ^0.05ML of Er is deposited, long 
chains of surface dimers form, followed by nanowires, 
with no intermediate reconstruction. Indeed, there are 
indications |120| that the surrounding substrate is clean 
Si(OOl). 

Another group of metals form closely related recon- 
structions: EuJi)8,J21j, YbW6, 118, lOSjUHl, NdfTUni and 
Sm |104| all form (2x3) and in some cases (2x4) recon- 
structions; by contrast to Er, Dy, Gd and Ho, these 
metals are strained in both surface directions on con- 
tact with Si(OOl) (which might be expected to reduce 
the periodicity of any substrate wetting) and lack the 
main driving force seen before for nanowire formation. 
In general these metals do not form nanowires [94111061 
I118j though there are certain conditions where Yb can 
be made to form nanowires |l()9j of some kind, which 
may well be grown off-axis (in a different growth mode 
to the other silicide nanowires which grow with the c-axis 
aligned with the substrate). 

It is also possible to induce nanowire formation in 
situations where they would not normally form, either 
because of uniform lattice matching (e.g on Si(lll) sur- 
faces, as mentioned above) or lattice mismatch in both 
surface directions using step edges. This technique has 
been successfully employed for Sm on SifOOLl ,90-92, and 
Gd[TTT)j and Dy[n2| on Si(lll). 

3. 2 Electronic Properties of Rare-Earth Nanowires 

One of the original technological interests in rare-earth 
overlayers on Si(lll) was the good crystal growth possi- 
ble due to the interface between the silicon and the over- 
layer, and the conductivity properties. The films have 
good conductivity and a low Schottky barrier (~0.4 eV 
on n-Si and ~0.8 eV on p-Si[HSl). 
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Measurement of the conductivity of RENWs is a sig- 
nificant challenge. Full two-point or four-point measure- 
ments requires either a unique STM instrument (for in- 
stance one used to measure conductivity of CoSi2 nanowires 
on Sifll0l [ll5| '). or alternatively the formation of nanoscale 
contacts to a NW |128j . While significant progress has 
been made in this area recently, any measurements of 
conductivity will include the resistance of the interface 
between the NW and the contacts. 

Scanning tunnehng spectroscopy (STS) has been per- 
formed on various NWs: Dy and Ho |MIITI7^ and GdjEl 
1^ . While STS does not measure the conductivity along 
the NW, it does measure the local electronic structure, 
and in all cases the NWs are found to be metallic, while 
the surrounding reconstructed substrate is not |94lll()8l 
I92| . This is good evidence that the NWs are taking 
on the bulk silicide structure, which is conducting, and 
the wetting layer is an intermediate state which relieves 
strain. 

Conductivity measurements of transition metal sili- 
cides on Si(lll) and Si(llO) have been made 1 12,115) . 
These studies exemplify the two techniques for measur- 
ing NW conductivity: NiSi2 NWs were contacted by gold 
padsinS, while four-probe STM measurements were made 
on CoSi2 NWs lllf^] . These NWs are relatively large com- 
pared to the RENWs discussed so far: 15nm (Ni) and 
60nm (Co) wide. The CoSi2 NWs showed a high Schot- 
tky barrier and conductivity equivalent to that of high 
quality thin films of CoSi2, while the NiSi2 NWs showed 
some signatures of quantum transport, though the con- 
ductivity of these wires was significantly lower than thin 
films; this is likely due to the overgrowth of the sample 
with Si02 for transfer to the lithography apparatus. 

3. 3 Reactivity of Rare-Earth Nanowires 

There is very little data on the chemical reactivity of the 
rare-earth nanowires. Thin films on Si(lll) are suscep- 
tible to reaction with 0[HH]; the RENWs oxidise rapidly 
in air jl2(Jj . and transition metal NWs show signs of this 
(confirmed by the conductivity of NiSi2 NWs on Sifllll fTT^ 
In the TM system, the NWs were overgrown with native 
oxide, and showed a significant decrease in conductivity 
relative to complete thin films of similar thickness, which 
is attributed to scattering at the NW surface/oxide in- 
terface. 

Further information on the reactivity comes from the 
deposition of Ft on a surface containing ErSi2 NWs on 
Si(001)'12U'. The Ft was deposited at room temperature 
and the sample was subsequently annealed. When STM 
images were taken they showed that the c(2x2) surface 
reconstruction on the NWs was no longer visible though 
the (2x1) substrate reconstruction was still present. Fur- 
thermore, the Pt-covered NWs were resistant to reactive 
ion etching and appeared stable when exposed to air for 
periods of up to 8 weeks; the Ft overlayer does appear to 



strain the NWs, however, leading to delamination of the 
NWs from the substrate if left untreated. Clearly there 
is much work to be done understanding the reactivity 
and stability of these NWs. 

4 Bi Nanovifires 

The discovery that nanowires would form if a Bi-covered 
Si(OOl) surface was left to anneal around the Bi desorp- 
tion temperature was made quite by chance. In 1995, 
in the Materials Department of Oxford University, two 
of the authors were investigating the surfactant-assisted 
growth of Ge/Si(001) using Bi as a surfactant, and were 
studying the properties of Bi/Si(001). When a sample 
which had been left to anneal at ca. 500°C overnight 
was imaged, very large, flat terraces, and long, straight, 
bright lines (longer than the scanning range of the STM) , 
running across the surface were found. These were the Bi 
nanolines. The first published image of the "nanobelts" 
was in 1997[S], and the first big study of their growth 
and properties came in 1999 9 . In this section, we will 
discuss the atomic structure of the Bi nanoline, and its 
physical and electronic properties. We show that many 
of the properties stem from its unusual structure. The 
reaction of the Bi dimers with a variety of reagents, and 
the effect of burial of the nanoline will also be discussed. 



Physical Structure 

The physical structure is the key to understanding the 
properties of a nanoline system, and great effort has been 
devoted to identifying the structure of the Bi nanoline. 
In this case, the structure was identified by a synergy 
of STM observations with tightbinding and DFT sim- 
ulations. The three different proposed models for the 
Bi nanoline are shown in Fig. |S1 The two early mod- 
els (shown as the top two models) share certain struc- 
tural motifs. They contain a pair of Bi dimers set into 
the surface layer of the Si(OOl) crystal, whose compres- 
sive stress is relieved by missing dimer defects (DVs). 
The first model proposed, shown in Fig.lSJa), was based 
around a pair of Bi dimers in the surface layer, sepa- 
rated by a missing dimer defect, so as to relieve the lo- 
cal stress of the Bi dimers pilll25| . This became known 
as the Miki model. The second proposed model, shown 
in Fig.E^b), was based around a pair of Bi dimers with 
defects on either side |126| . which became known as the 
Naitoh model. However, the third model is more com- 
plex: a reconstruction of several layers of Si underneath 
the Bi dimers produces the nanoline core. This was named 
the Haiku model [T77| ^. Of the three models, this struc- 
ture is the only one to fit all the criteria which have been 

It is named after the Japanese verse form, which com- 
prises three lines of 5, 7, and 5 syllables, referring to the 5- 
and 7-membered rings of silicon in the reconstruction. 
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Fig. 5 Ball-and-stick models of the three proposed models 
for the Bi nanoline. (a) The Miki model. 0. (b) The Naitoh 
modelfn?)]. (c) The Haiku model[T?f]. 



determined from STM observations, as well as having 
the lowest energy of the three. In this section, we will 
briefly review the historical process by which this struc- 
ture was determined, and discuss the properties of the 
Haiku structure. 

J^.1.1 Early Models Early STM images of the Bi nano- 
line were taken at elevated temperature pill25| . in which 
the Si dimers were not resolved, and hence the registry 
of the nanoline with the substrate was unknown. This 
data suggested that the nanoline contained two features, 
which were probably Bi dimers, with a spacing of ap- 
proximately 6.3A, with a total width of about 1 nm, 
which was approximately equivalent to the space of three 
Si dimers. As the nanolines appeared bright at large bias 
voltages (in both positive and negative bias images) to 
dark at low voltages, it was also reasonable to assume 
that these Bi dimers were situated in the surface layer, 
rather than in an adlayer or in a subsurface layer. It was 
thought that Bi dimers embedded in the top surface layer 
would have considerable compressive stress, and on this 
basis, a model based around a IDV, with Bi dimers to ei- 
ther side, was proposed. This was the Miki modei pill25[ 
inH]. Tightbinding calculations of the Miki modelfHSI 
found that it was more stable than the (2x1) or (2xn) Bi 
reconstructions, and that the formation energy of a de- 
fect in the line was high: around 1.1 eV |128| (this number 
considers putting the Bi which has been removed as an 
ad-dimer on the surface; recent data |129] shows that the 
defect energy falls to 0.11 eV if the Bi dimer is placed in 
another Miki model; this point is discussed more fully in 
Section ETl below). DFT calculationS tl30a27...35 ,131^ 
agreed with those conclusions. The calculated LDOS 
showed that the Bi dimer states were further away from 
the Fermi level than those of the Si dimers and simulated 




Fig. 6 10 nmxlO nm STM image of the Si(001):H surface. 
The H termination makes it easier to resolve the individual 
Si dimers. A series of markers has been placed across the Bi 
nanoline. From this, it can be determined that the nanoline 
occupies the space of four dimers. 



STM images of the Miki structure [ 



showed 



that the line would appear dark at low bias voltages, in 
agreement with STM results [T5^133| . Testing by other 
methods appeared to confirm this structure. Photoemis- 
sion spectroscopy experiments T27' found that the Bi 5d 
core- level spectra of the Bi nanowire was essentially iden- 
tical to the spectra of the (2xn) phase composed of Bi 
ad-dimers. This suggested that the local chemical state 
and registry of Bi adsorbates for both phases was the 
same, i.e. that the Bi was in the form of dimers in the 
top layer of the structure. X-ray photo-electron diffrac- 
tion (XPD) experiments |134| found a good fit between 
the experimental XPD data and simulated XPD inten- 
sity peaks from the Miki model. In particular, they con- 
firmed the presence of Bi dimers parallel to the Si dimers 
and found the spacing between them to be '-~^6.3A, in 
agreement with the STM measurement, though it is im- 
portant to note that the quoted distances are the result 
of a fit to the Miki model. 

However, it was difficult to account for the extreme 
straightness of the nanoline with the Miki model, as the 
calculated kinking energy was very small, around 0.1 eV. 
A kinetic argument was put forward, which suggested 
that there was a stronger preference for incoming Bi 
dimers to line up at the end of the nanoline, and hence 
the nanoline would grow straight j^. Furthermore, the 
diffusion constant of a Bi surface dimer was high (this 
diffusion was required for a kink to form in a Miki model 
nanoline), and hence a line would remain straight once 
it had grown. Despite this explanation, the issue of the 
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straightness of the Bi nanolines remained as a question 
mark over the Miki model. 

Subsequent room-temperature STM experiments j^Bj 
proved decisively that the Miki model could not be the 
structure of the Bi nanoline, as it had the wrong registry 
with the surface. Markers were placed on the surface of 
an image with the Si dimers resolved, as shown in Fig.O 
and by counting across the nanoline, the width was de- 
termined definitively to be four dimers, and not three as 
in the Miki model; moreover the Bi dimers lay between 
the Si dimers of the substrate, not in the same posi- 
tion. (This registry has been confirmed by every other 

On the basis of the new STM data, the Naitoh model was 
proposed |I126I| . Again this model was tested by atomistic 
calculations; a detailed study of the dimensions and sim- 
ulated STM appearance llHlUflfiUlHS ' showed that again 
the nanoline would appear dark in STM at low bias volt- 
ages, but showed also that the spacing of the Bi dimers 
was too narrow, approximately 5A, which could not be 
reconciled with the spacing of 6.3A measured from STM. 
Moreover, the energy was considerably worse than the 
Miki model [TTTItnT] . and the kinking energy was still 
low. There was therefore no satisfactory structural model 
at this time. 

4-1-2 The Haiku model Inspiration for the structural 
model of the Bi nanoline arose from the structural model 
proposed for the B-type double-height step of As-terminated 
Ge |139| . At this step edge, bond rotation forms a pair of 
5 and 7-membered rings of Ge, capped with As, shown 
in Fig. Ufa) . While the Haiku structure appears to be a 
complex reconstruction, in fact the construction from a 
cell containing two Bi ad-dimers can be described schemat- 
ically in two simple steps, as shown in Fig.|3b)-(d). The 
first step is to rotate the second and third layer atoms 
underneath the Bi ad-dimers, so that the Si dimer atoms 
sink down into the surface. This creates the 5- and 7- 
membered rings also seen in the As/Ge step structure. 
The central four atoms are then removed, and the two 
halves of the structure bonded together, thus creating 
the Haiku structure. (N.B. This is not proposed as the 
formation mechanism, but is simply given as a means 
of understanding the substructure.) A linescan of the 
nanoline has been matched to the Haiku model, using 
the background Si dimers as reference marks [T?7| . The 
positions of the Bi dimers match up extremely well to 
the peaks in the linescan, giving strong confirmation 
of this model. Other features of experimental linescans 
also agree well with simulated linescans for the Haiku 
model [nnmSZI- DFT calculations of this new structure [T?7] 
found that it was considerably more stable than the 
Miki model, which has been confirmed by subsequent 
modelling |14lll33lll37lll38j . Although the Haiku struc- 
ture shares several structural features with previous models— 
the top of the structure has two Bi dimers joined by 
rebonded second-layer Si atoms, as in the Miki model, 




Fig. 7 Ball-and-stick models of the As/Ge step edge struc- 
ture, a pair of ad-dimers, the double-core structure, and the 
Haiku structure. Schematically, a double-core structure may 
be achieved by rotation of the atoms marked 'A' underneath 
the ad-dimers from a vertical plane into a horizontal plane. 
Removal of the 4 central atoms marked 'B' and rebonding 
then produces the Haiku structure. The 5-7-5 structural mo- 
tif is also marked. 

while on the outside of the Bi dimers, there are more 
rebonded Si atoms, just as in the Naitoh model — the 
difference lies in the Si substructure; this mixture of 5 
and 7-membered rings extends down many layers. Con- 
sideration of the core of the nanoline reveals that it is 
a small triangular section of hexagonal silicon, embed- 
ded in the diamond cubic silicon substrate, rather like 
the "endotaxial" RE nanowires mentioned in Section l3.ll 
The {111} planes that delineate the core are shown in 
red in Fig. |H1 The hexagonal core of the nanoline does 
not exist as a bulk structure for silicon; it is particu- 
lar to this nanoline structure. In this respect, the Bi 
nanoline is quite unlike the epitaxial rare-earth silicide 
family, which are essentially one bulk crystal grown epi- 
taxially onto another bulk crystal, taking their shape 
-from the highly anisotropic mismatch between the two. 
This core structure is responsible for many of the prop- 
erties of the Bi nanoline. It will be difficult or impossi- 
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Fig. 8 A cross-sectional perspective view of the Bi nanoline. 
The nanoline is 1.5 nm wide, occupying the space of 4 surface 
dimers, and is built around a pair of Bi dimers, on top of a 
complex Si substructure. The shaded {111} planes define the 
interface between the nanoline core and the silicon substrate. 



ble to kink, resulting in the extremely straight nanoline 
observed. As will be shown below in Section [4.2.41 the 
tensile strain field results in the repulsion of other nano- 
lines, step edges and missing dimer defects. As a result, 
the nanolines will not grow sideways, will not grow to- 
gether, and will not coarsen into larger islands, unlike 
the silicide wires. The electronic structure of the Haiku 
model has been studied in some detail [TTi2ll 1 40111 HHU 1 88| . 
and again the nanoline becomes dark in low-bias images 
— this effect is a property of the Bi dimers rather than 
that of the nanoline. Low-bias STM images reveal other 
electronic effects, which provide further confirmation of 
the Haiku structure, while ruling out the Miki structure. 
These features are discussed below in Section ^31 

However, despite the wealth of agreement of the pre- 
dictions from the Haiku model with experimental evi- 
dence, there has been no direct confirmation of the Haiku 
structure. The Bi nanolines are destroyed by exposure to 
air, while burial of the nanolines in a capping layer also 
changes their atomic structure |141lll42| (though if done 
carefully the ID character and dimerisation of the Bi can 
be preserved), so that cross-sectional TEM studies [T^ 
have not provided definitive confirmation of the core 
structure, and a recent X-ray Standing Wave(XSW) study 
which relied upon capping with amorphous silicon and 
whose results cast doubt on the Haiku model, cannot be 
regarded as a counterindication. The difficulties of burial 
of the Bi nanoline are discussed in Section IT^ We note 
also that a recent suggestion |188j . that the height of the 
Bi dimers in the Haiku model is too far above the plane 
of the surface Si dimers based on XPDjIJJ, is not a sig- 
nificant piece of evidence against the Haiku model: the 
XPD data was fitted to the Miki model, thus naturally 
giving the wrong data for the Haiku. If a fit of the data 
to the Haiku model were to give height data which con- 
tradicted the atomistic modelling, then the model would 
have to be revisited. 



^. 1.3 Miki model revisited Despite the apparent suc- 
cess of the Haiku model, the Miki model has contin- 
ued to attract a large amount of attention 130 131 , 5^ 
I137II138| as a possible candidate for the nanoline struc- 
ture. While in fact its incorrect registry with the Si sub- 
strate unambiguously rules it out as a candidate nanoline 
structure |12f)lll27lll29| . it is an energetically favourable 
structure, and might be expected to be present on the 
Bi-covered Si(OOl) surface. With the benefit of hindsight, 
it can be seen that the Miki structure was indeed ob- 
served in early STM images of the Bi-rich surface. Fig.|51 
shows that annealing a Bi-rich surface at 400° C pro- 
duces a large density of quite straight dimer vacancy 
(DV) trenches. In the filled-states image, Fig.|5Ja), these 
appear to be composed of clean silicon. However, in the 
empty-states image, Fig.l^Jb), there is contrast between 
the dimers adjacent to the missing dimer trenches and 
the rest of the surface. This contrast docs not occur on 
the clean surface at this bias voltage, and so it can be 
inferred that the trenches are decorated by Bi dimers, 
i.e. the Miki structure has been formed. As with the Ge- 
induced DV trenches mentioned in Sec.|3| there is coop- 
erative strain relief between a Bi dimer and a IDV. Bi 
dimer/DV structures have the lowest energy/Bi dimer of 
any structures except the Bi nanolines. There is a kink- 
ing energy for these (2xn) trenches of around O.leV fT^ . 
Hence the formation of a semi-regular array of Bi dimers 
reduces the energy further. 



Further proof of the Bi decoration of IDV trenches 
has been provided by a more recent experiment, in which 
a Bi-rich surface was quenched from 550 °C at an early 
stage of the annealing process 129, . Here atomic hydro- 
gen was adsorbed after cooling to enhance the contrast 
between Bi and Si significantly, as the H will adsorb 
easily to the Si dimers, but does not adsorb on Bijl2^ 
I135II129| . as can be seen in Fig. |H1 It was found that 
the H termination increased the relative apparent height 
if Bi-related features by ca.l50 pm, or 1.5A. An ex- 
mple image is shown in Fig. Efc). A linescan across 
both a Bi-decorated trench and a section of nanoline 
shows that the Bi-decorated trench is quite similar to 
the nanoline, but is lower and has a deeper depression 
between the two Bi dimers. STM linescans of the fea- 
ture decorating DV trenches and the Bi nanolines |129| 
compare very well to simulated scans of the Miki model 
and Haiku model respectively^^, providing strong sup- 
port for both these identifications. Furthermore, in this 
way, it has been shown that the Miki structure and the 
Bi nanoline co-exist, demonstrating once again that the 
Miki structure is not the nanoline, as has been recently 
suggested p^s, but it is kinetically stable, even at high 
temperatures, due to the large activation barrier neces- 
sary to form the Haiku structure. 
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Fig. 9 STM images of Si(OOl) surface after deposition 
of Bi for 70s at 400 °C, followed by 10 min. annealing, 
(a): -2.2 V,0.3 nA. There are many straight missing dimer 
trenches, (b): +2 V, 0.3 nA. Brightening around the trenches 
in empty-states images indicates that they are decorated by 
Bi dimers, i.e. the Miki structure, (c): A Bi-decorated trench 
and Bi nanoline on a H-terminated surface. A linescan along 
the dotted line (d) reveals that the Bi-decorated trenches are 
significantly lower than the nanoline. (a) & (b) are reprinted 
from Appl. Surf. Sci. 142, Naitoh et al. Bismuth-induced sur- 
face structure of Si (100) studied by scanning tunneling mi- 
croscopy, p. 38, (1999), with permission from Elsevier. |145| 



4-. 2 Formation and Properties of Bi nanolines 

While rare-earth silicide wires have quite specific growth 
recipes in order that the nanowires do not coarsen, the 
Bi nanolines are robust to a wide variation in the growth 
conditions. The essence of Bi nanoline formation is that 
it is a competition between deposition and evaporation. 
A thick layer of Bi on Si(OOl) is only stable below the Bi 
bulk evaporation temperature (ca. 400°C). Above this 
temperature, an epitaxial layer of 1-2 ML Bi is stable, 
and annealing of this surface will result in a surface with 
Bi nanolines co-existing with these Bi islands, as shown 
in Fig ^] The first-layer Bi forms a (2xn) structure, 
with n=4 or S pPUT^ . and the second-layer Bi forms as 
small groups of dimers, which can be either parallel or 
perpendicular to the underlying Bi. 

The threshold temperature at which Bi dimers on 
Si(OOl) will start to evaporate is around 500°C, while 
the maximum temperature at which any Bi is stable 
on the Si(OOl) surface is around 600°C. Deposition of 
Bi in this temperature window, or deposition of Bi at 
a lower temperature, followed by an anneal within this 
temperature window, will result in a surface comprising 
Bi nanolines on an otherwise clean surface, as in Fig. II II 




Fig. 10 STM images of silicon surface onto which 3 ML of 
Bi has been deposited at 360° C, and then annealed for 60 
mins at 497° C, resulting in Bi nanolines surrounded by a Bi 
overlayer, which forms a (2xn) reconstruction, with n=4 or 5, 
to accomodate its large strain. Left-hand images are empty- 
states, right-hand images are filled-states. The (2xn):Bi does 
not grow over the Bi nanoline at this temperature. 



However, the details of the growth recipe will determine 
the surface morphology at the end of the anneal, as is de- 
scribed below in Sec. 14.2.11 Of particular technological 
interest is the growth of majority-domain surfaces af- 
ter long anneals. The Si(OOl) surface has two equivalent 
domains with the dimer rows running in orthogonal di- 
rections. On a typical Si(OOl) surface, these two domains 
will have roughly equivalent total areas. Single-domain 
surfaces form on vicinal wafers miscut along a (110) di- 
rection, where the dimer rows tend to run perpendicular 
to the step edges. However, by long anneals of the Bi- 
rich surface towards the high end of the temperature 
window, majority-domain surfaces can be formed, even 
on flat surfaces. In this case, extremely long Bi nanolines 
grow, up to 1/xm, as can be seen in Fig. 1121 

High-temperature STM observations of the annealing 
Bi-rich surface shows a consistent series of surface mor- 
phologies. Immediately after the end of Bi deposition, a 
relatively rough surface is obtained, with a high density 
of small islands, composed of a mixture of Bi and Si, as 
shown in Fig. Illf a). These will quickly disappear, but 
at the same time, nanolines begin to form, either on the 
terraces, or on one of the mixed islands. The nanolines 
grow out over the step edges in long, finger-like islands 
known as promontories, with a Bi nanoline surrounded 
by a thin strip of silicon. Examples of the growth of 
promontories are seen in Fig. 1171 and more mature ex- 
amples of this type of surface are shown in Figs.lTTiyT^ 
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Fig. 11 Different growth procedures for the formation of Bi 
nanohnes, and the resulting surfaces, (a): Deposition of 1 ML 
Bi at a lower temperature, followed by a 20 min. anneal at 
520° C. (b): Deposition of 1 ML Bi at around 550° C, followed 
by annealing at the same temperature, (c): Continuous ex- 
posure of 0.1 ML/min Bi flux for 40 mins at 570°C. (d): As 
(c), but on a 0.5° miscut sample. The growth recipes are 
summarised schematically in Fig |13l 

The nanolines will not only form long, narrow promonto- 
ries, but also the reverse, a deep inlet where a nanoline in 
one terrance has grown through an up-step into an upper 
terrace. The surface between the Bi nanolines contains a 
large density of mobile missing-dimer defects, apparent 
in STM as dark, fast-moving streaks. After further an- 
nealing, the missing dimer trenches begin to disappear, 
and the Bi nanolines remain on a very flat, featureless 
surface. Finally, the Bi nanolines themselves evaporate. 
In the following sections, we will give describe the sur- 
face processes which drive this sequence of events, and 
account for these observations. The interactions of the 
Bi nanoline and its associated strain field with defects 
and step edges in the terraces, and with each other, will 
be discussed and the possible nucleation mechanisms of 
this complex structure explored. 

4-2.1 Growth Recipes As shown above, annealing a Bi- 
rich surface at lower temperatures results in a surface 
containing a high density of missing dimer trenches, dec- 
orated with Bi dimers, i.e. the Miki structure. One fea- 
ture of the Bi nanoline is the high temperature required 
to form it, implying a significant activation barrier to its 
nucleation. While the large number of growth param- 
eters, such as flux, total amount deposited, substrate 
temperature, miscut angle and azimuth, provides a wide 
range of different possible recipes, as summarised in Fig. 1131 
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Fig. 12 (a): An example of a single-domain surface with Bi 
nanolines hundreds of nm in length. Very long promonoto- 
ries and inlets are present. Reprinted with permission from 
Ref Copyright (2004) by the American Physical So- 

ciety, (b): The larger-scale image shows clearly the single- 
domain nature of the surface, with double-height steps sep- 
arating domains. Images courtesy of Prof. A.B. McLean, 
Queen's University, Kingston, Ontario, Canada. 

there are essentially two routes to the growth of nanolines^. 
One growth method for the nanolines is to deposit a 
monolayer of Bi below this threshold temperature, and 
then anneal the sample until Bi nanolines form. After 
around 1 hour at around 500° C, this will result in a 
surface with Bi nanolines and a background Bi-covered 
surface |8u9,. Continued annealing of this surface will re- 
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suit in the desorption of the background Bi forming a 
surface similar to that shown in Fig lllf a'). and even- 
tually a surface similar to that shown in Fig. Ilir b) is 
obtained. The second growth preparation route involves 
the deposition of Bi onto the surface within the des- 
orption temperature window. Due to the lower surface 
coverage of Bi, and the higher temperature, a more or- 
dered surface is obtained, and fewer, longer, Bi nanolines 
form. The final surface is much the same, however, cor- 
responding to that shown in Fig llir b'l. A high density of 
Bi nanolines can be built up by continuous Bi exposure 
for a long time, at a temperature in the upper half of the 
desorption window. In this regime, the nanolines are only 
stable for a few minutes in the absence of a Bi flux and 
small differences in stability, such as the ends of a nano- 
line, will be significant. Shorter nanolines, perpendicular 
to the prevailing direction, will tend to dissolve, so that 
a highly ordered surface with long, parallel nanolines is 
obtained, as in Fig lllf c^. By use of a vicinal surface, 
the preference for one domain can be enhanced further, 
and a surface close to a single-domain surface can be ob- 
tained, as in Fig llir dl. A second method which produces 
a majority-domain surface, even on a surface without 
any intentional miscut, (the miscut angle is (0.5°) 
I146| . produces a surface as shown in Fig. ^1 In this 
case, a substrate temperature around 530° C was used, 
4.5 ML of Bi were deposited, with an anneal for around 
40 mins. Comparison of Fig lllf d^ and Fig. El reveals 
that in the former case, the nanolines are parallel to the 
prevailing step edges, while in the latter, they are per- 
pendicular. Despite the large difference in morphology, 
the only differences between this recipe and that which 
produced Fig lllf c.d'l are that continuous deposition for 
an extended period was not used, and the heating cur- 
rent used was AC, a strategy which was designed to re- 
duce electromigration. Whether this is is the cause of 
the majority-domain surface remains unclear, but such 
large-scale single-domain surfaces are ideal for nanoelcc- 
tronics applications. 

While STM can observe only a relatively small part 
of the surface, the large-scale order of the surface can 
be determined using electron diffraction methods such 
as RHEED. A series of RHEED patterns at different 
azimuthal angles of a single-domain Bi nanoline surface 
are shown in Fig. 1141 Perpendicular to the nanolines, the 
RHEED pattern comprises vertical lines, while parallel 
to the nanolines, a ring feature is seen. The strength of 
these features demonstrates the long-range order gener- 
ated by the nanolines, and means that their growth and 
development can be followed in real-time. 

4-2.2 The Bi reservoir In contrast to the silicide nanowires, 
which form very rapidly at the growth temperature and 
are only stable for about 10 mins. at that temperature, 
there is a significant incubation period before the nucle- 
ation of the Bi nanolines, and they will continue to grow 
for upwards of an hour after the deposition of Bi has 




Fig. 13 Schematic diagram showing the different recipes 
used to produce Bi nanolines. Below 500° C, a ML of Bi is 
stable on the surface, and the (2x5):Bi reconstruction re- 
sults. Annealing this surface above 500° C, or deposition of 
Bi onto the clean surface above this threshold results in Bi 
nanolines on a clean background. High coverages of nanolines 
are obtained by higher-temperature anneals, with longer Bi 
exposure times. 



finished. Thus the Bi required to generate the nanolines 
must come from a reservoir of Bi which is present on the 
surface. Where is this Bi located, and how does it affect 
the growth of the nanolines? 

In images of the surface taken soon after the end of 
the Bi deposition, for example, Fig llir aL there is often a 
high density of small islands on the surface. These islands 
have the same appearance as the rest of the surface; they 
are not islands of Bi dimers, as may be seen after deposi- 
tion at lower temperatures, but appear to be mixtures of 
Si and Bi atoms. In general, such small islands are not 
stable at these temperatures, as has been seen in hot 
STM studies of Si homoepitaxy on SifOOP pil^ . Dur- 
ing growth, these small islands are stabilised by the high 
flux of Si atoms, but if the flux is cut off, a process similar 
to Ostwald Ripening occurs, in which these islands will 
decompose, and the material in them will move to step 
edges (which may be regarded as a very large island). 
The similarity of the behaviour observed here suggests 
that these islands are sustained by a flux of Si and Bi, 
which is produced during Bi deposition. For low cover- 
ages of Bi, the lowest-energy sites for Bi dimers which 
have substituted for Si dimers in the top layer of the 
Si(OOl) surface — hereafter "Bi surface dimers" — are Bi 
surface dimer/missing-Si-dimer complexes |1 29j . where a 
Bi dimer decorates one or both sides of a IDV, as seen 
in Fig. 1^1 For each adsorbed Bi dimer, therefore, up to 
four Si atoms are ejected from the surface layer, pro- 
ducing a significant transient flux of Si atoms; this flux 
in turn produces the observed islands. At the end of Bi 
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Fig. 14 A series of RHEED images of a Bi nanoline surface. 
Along the [110] direction, perpendicular to the nanolines, a 
series of long streaks is seen, as in (a). Along the [110] di- 
rection, parallel to the nanoline, as in (c), a distinctive ring 
feature is seen, while at intermediate angles (b), the ring 
opens up into a long arc. 



deposition, this flux will die away, the islands will no 
longer be stable and will decompose, as occurs between 
Fig llTf aland (b). A cartoon of this process is shown in 
Fig. cad). 

Once the small islands have annealed away, the sur- 
face morphology is very flat, and the background be- 
tween the nucleating Bi nanolines appears featureless in 
hot STM, apart from missing dimer trenches. However, 
at this stage a large surface density of Bi may be inferred 
from the fact that Bi nanolines continue to nucleate and 
grow on this surface. The Bi dimers embedded in the 
surface layer have approximately the same contrast as 
the Si dimers [T^ . and are therefore invisible at elevated 
temperature, as in Fig. I15f a). At a large positive bias, 
however, the background Bi becomes visible, as is shown 
in Fig. llSf b). Top-centre of this image is a Bi nanoline. 




Fig. 15 (a): Hot STM image of the mature nanoline surface, 
showing Bi nanolines and dark missing dimer defect rows. 
Any background Bi is invisible, (b): Hot STM image of the 
Bi:Si(001) surface taken at about 500° C, with a sample bias 
ca. +3.0 V. Grey linear features, which run parallel to the 
Bi nanolines, can be seen in the background. An example 
is marked by white arrows, (c): STM image taken at room 
temperature, with the silicon passivated with hydrogen to 
increase the contrast between the Si and Bi. The grey linear 
feature can now be seen to comprise Bi dimers decorating a 
missing-dimer trench, i.e. the Miki structure, (d): A cartoon 
of the formation of the Bi reservoir. Impinging Bi atoms (1) 
exchange with Si atoms, making Bi surface dimers and Bi!DV 
complexes (2). The ejected Si atoms difffuse and form islands 
(3). 



There is a high density of linear grey features in the 
background, one of which is marked with a pair of white 
arrows. Detail of these features cannot be discerned at 
this temperature, although the background Bi is notice- 
ably lower than the nanoline. After quenching the sample 
to room temperature, and exposure to hydrogen, we can 
judge the true coverage and distribution of Bi during 
the annealing process, as displayed in Fig. llSf c). There 
is a high density of Bi surface dimers, and also many Bi 
dimers are adsorbed to one or both sides of a missing 
dimer trench, forming the BilDV or Miki structures, as 
shown in Fig. |5| Thus the Bi reservoir consists not of 
mobile ad-dimers, but mostly of Bi surface dimers, usu- 
ally decorating the missing-dimer trenches to form the 
Miki structure. 

During the high-temperature anneal, the surface Bi 
reservoir is continuously depleted, with some Bi going to 
form nanolines, and the remainder evaporating back into 
vacuum. The activation barrier for evaporation is dif- 
ferent for different surface adsorption sites. A RHEED 
study[2 found that the activation barrier for evapora- 
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Fig. 16 A sequence of hot STM images, showing the 
breakup and evaporation of the BI nanolines at high tem- 
peratures. Arrows indicate the major areas of change from 
one image to the next. 



tion of Bi from the nanoline was 0.25 eV higher than the 
barrier for evaporation of Bi from a (2xn):Bi overlayer. 
Hence Bi is mostly being lost from the background reser- 
voir. However, it might be expected that missing dimers 
in the nanolines would be observed occasionally. In fact, 
this is not the case, and those nanolines which are ob- 
served at high temperatures always appear perfect. It is 
likely that, due to the large defect energy of the Haiku 
structure (0.66 eV from DFT calculations T5I ) , any Bi 
which evaporates from the nanoline is quickly replaced 
from the Bi reservoir. Thus at high temperatures, where 
Bi can evaporate from the nanolines, their stability is 
dependent upon a sufficient quantity of Bi in the reser- 
voir to replace missing dimers. At the later stage of an 
anneal, no new Bi nanolines nucleate, and those which 
are present cease to grow. This is an indication that the 
reservoir of Bi is approaching exhaustion. Finally, when 
the reservoir is exhausted, evaporation from the Bi nano- 
lines will take place, and they will disappear within a few 
scans, or sometimes within the time of one scan. A se- 
quence showing the evaporation of nanolines is shown 
in Fig ll6l In general, the nanolines evaporate from the 
ends, rather than breaking up into pieces, which strongly 
suggests that the ends of nanolines are less stable than 
the middle of the nanoline, and that the energetic bar- 
rier to break the nanoline in the middle and create two 
ends is prohibitive, even at these temperatures. 

4-2.S Promontory Growth The most notable feature of 
the mature nanoline surface as seen in Fig. Illf c) and 
(d) as well as Fig. ^1 is known as the promontory. Long, 
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Fig. 17 A sequence of high-temperature STM images show- 
ing the growth of a promontory. Between the first and second 
images, the high density of islands has annealed away, and 
two nanolines have nucleated. The growth of nanolines in 
the succeeding image is indicated by the white arrows. The 
promontory is a growth efi'ect, and not caused by etching 
away of the surrounding Si. 

narrow islands of Si surrounding Bi nanolines grow out 
across lower terraces. A series of images showing the 
early stages of Bi nanoline growth are shown in Fig. 1171 
In (a), a short nanoline has grown out over a lower ter- 
race and the joining of the promontory and the pre- 
existing step edge has formed a hole in the upper terrace, 
as is marked in (b). The growth of this promontory is 
blocked by the growth of a second nanoline, which has 
nucleated in the time interval between (a) and (b), and 
is growing out across the next lower terrace. Compar- 
ison of (b),(c) and (d) makes it it very clear that the 
promontories are the result of growth out over a lower 
step edge. By (d), a common pattern of zigzagging Bi 
nanolines which have blocked each other's growth is de- 
veloping. 

The mechanism behind this promontory growth has 
not been determined; in an early paper, it was suggested 
that the promontories form by etching of Si around the 
Bi nanolines However, in Fig. El there is no sign of 
etching taking place; this is a pure growth phenomenon. 
It is likely that the growth of the nanoline is related 
to the high mobility of step edges at 500° C. Nanolines 
will always continue to grow until they reach an obsta- 
cle. Thus nanolines will stop at step edges. However, 
the position of the step edge at this temperature is not 
constant; it fluctuates back and forth |14'7) . as material 
moves randomly along it. This motion may be inferred 
from the jagged, streaky step edges in these high-temperature 
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Fig. 18 Left-hand side.: STM sliowing parallel Bi nano- 
lines. The positions of missing dimer defects in between the 
two nanoline marked with white arrows have been gathered 
and are shown ini the graph below. (Reprinted from Appl. 
Surf. Sci. 412, Naitoh et al. Bismuth-induced surface struc- 
ture of Si (100) studied by scanning tunneling microscopy, 
p. 38, (1999)123, with permission from Elsevier.) Right-hand 
side: Graph of interactions between nanoline and other sur- 
face features - IDV defects, step edges, Bi surface dimers, 
and other nanolines. Distances are measured from the centre 
of the nanoline to the centre of the feature. High-temperature 
image showing Bi nanolines and disordered dark trenches of 
Si missing dimers. Although the defects are very mobile, there 
is an empty space around each nanoline, known as the Defect 
Exclusion Zone, or 'DEZ'. The DEZ is about 7 nm across. 
This width is the same as the width of the promontory. 



images, as discussed in Section ITTl Therefore, each time 
the step edge moves forwards, the nanoline can extend 
to the new step position, and then when the step edge 
tries to move back again, its position is pinned by the 
presence of the nanoline. With each forward fluctuation 
therefore, the nanoline will grow forward, and thus a 
promontory will gradually form. A similar process could 
explain the retreat of an up step ahead of a growing 
nanoline. The notable feature of the promontory is that 
they have a fixed minimum width. At least 3-4 nm of 
Si is maintained on either side of the Bi nanoline. This 
fixed width is an example of the phenomenon known as 
the Defect Exclusion Zone, or "DEZ", which is discussed 
in Section 

4-. 2.4 The role of strain in surface interactions Since 
the nanoline has a strain field associated with it, it is 
likely to exhibit quite strong interactions with other sur- 
face features which exert a localised strain on the sur- 
face. Such features might be adsorbed Bi dimers, missing 
dimer defects, step edges, and other nanolines. 



Using tight-binding simulations, where very large cells 
can be relaxed, the interaction between a nanoline and 
several features, such as IDV defects and rebonded B- 
type step edges have been calculated. These data are 
plotted in Fig^J From elevated-temperature STM im- 
ages, it is clear that there is a repulsive interaction be- 
tween the nanolines and missing dimer defects and step 
edges. At high temperature, as in Fig. ^1 the surface 
is covered in black streaks, which indicate the positions 
of rapidly moving missing dimer defects. However, ei- 
ther side of the nanoline, there are no streaks. Likewise, 
where a nanoline forms a promontory, the width of this 
empty area is continued as the width of the promontory. 
This characteristic width of defect-free silicon, around 3- 
4 nm either side of the nanoline, is known as the "Defect 
Exclusion Zone" or DEZ |140| . For a large separation, 
the total energy of a nanoline and a defect is the same 
as for the two structures independently. However, as the 
separation decreases, the total energy increases, so that 
at a distance of about 3.5 nm, there is a repulsive inter- 
action between the nanolines of ca. 0.1 eV for both step 
edges and IDVs, and at smaller distances, this energy in- 
creases rapidly. Statistics about the position of IDVs in 
the space between two nanolines have been gathered, as 
shown in Fig. ^1 taken from Ref. |145| . They show that 
on average, most defects occur at least 3.2 nm away from 
the nanoline, which is in agreement with the calculated 
repulsion. 

For a pair of nanolines, there is a similar interaction, 
and the 0.1 eV threshold is also reached when the centres 
of the two nanolines are about 3.5 nm apart, or in other 
words, the gap between the nanolines is the width of one 
nanoline. Clusters of nanolines with this approximate 
spacing are often seen, as in Fig.^| For nanolines adja- 
cent to each other, the excess energy is about 1 eV per 
unit cell. Despite this, nanolines are sometimes observed 
growing next to each other, as in Fig. lllf c). The reason 
for this is that unlike defects and step edges, which are 
mobile, a nanoline cannot move sideways. Hence when 
two nanolines which have nucleated in different places 
grow past each other, they cannot move sideways to re- 
duce their interaction energy. The presence of nanolines 
in close proximity does indicate that this interaction en- 
ergy is not so high as to present an insuperable barrier 
to the growth of two nanolines past each other. For Bi 
surface dimers, which have a compressive stress field, 
decoration of missing dimer trenches is a natural way 
in which to relieve their stress. The gain in energy by 
this process is about 0.45 eV/Bi dimer in tightbinding 
calculations ^S]. The interaction with the nanoline is 
also attractive. In images where there is a large surface 
density of Bi dimers, they are often situated adjacent 
to the nanolines. This position is found to be 0.28 eV 
better in energv ,129) than elsewhere. This suggests that 
the nanolines are generally more stable in a surface with 
a significant amount of compressive stress, such as one 
which is covered in Bi. Indeed the energy of a nanoline 
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Fig. 19 Example of observation of the nucleation of the Bi 
line at 497°C. A series of 81 nm STM images. During the 
interval from (b) to (c) of 14 min., the Bi line on the bot- 
tom extended to separate the big Bi island and to make an 
inlet as marked by a white arrow. During the same period, a 
new Bi line appeared as marked by a black arrow, (d) to (e) 
shows more details of observation corresponding this latter 
event, in the area marked by a white box in (c). The time 
between images was 18 sec. In figure (e) a small Bi island 
is situated at the place marked by a black arrow, in (f) the 
small Bi island changed into a short Bi line, and between 
(f) and (g) the new line extends in both directions towards 
the island and the inlet. It is noted that the side edges of 
the inlet, and the shapes of all the islands change image by 
image. Reprinted from Surf. Sci. 421, Miki et al. Bi-induced 
structures on Si(OOl), p. 397, (1999), with permission from 
Elsevier. 



in a Bi-terminated surface is 0.8 eV better than on a 
Si-terminated surface. 

4-2.5 Nucleation Despite the success of the Haiku struc- 
ture in explaining all of the above surface phenomena, 
its identification raises as many questions as it answers. 
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Fig. 20 A proposed mechanism for the nucleation of a Bi 
nanoline. A Bi ad-dimer fills a BilDV trench, to form a pair 
of adjacent Bi dimers. In order to relieve its strain, the atoms 
between the Bi dimers twist sideways to form a zigzag struc- 
ture. Significant reconstruction takes place between (iv) and 
(v), resulting in the Haiku structure. The barrier for this 
process is less than 1.5 eV. Reprinted with permission from 
Ref. |15| . Copyright (2005) by the American Physical Society. 



How would such a complex structure form? How is it 
terminated? If termination is energetically unfavourable, 
what would a nucleus structure look like? 

Experimentally, the atomistic nucleation process re- 
mains a mystery. There is an incubation time during 
the annealing process before nanolines start to nucleate, 
which suggests that there is a significant barrier to for- 
mation of these structures. On many occasions, the ap- 
pearance of Bi nanolines has been captured during STM 
observations, such as in Fig. ^| There is no apparent 
precursor state; like Athena from Zeus' brow, a short 
nanoline springs fully-formed from an empty patch of 
the surface. The nanolines then grow extremely rapidly. 

Thus far, there has only been one proposed mecha- 
nism for the nucleation of a Bi nanoline|15|. It has been 
suggested that the initial step is for Bi ad-dimers to fall 
into the missing dimer trenches on the hot Si(OOl) sur- 



20 



J.H. G.Owen et al. 




face, forming a row of pairs of Bi dimers. Such a structure 
would result in the increase in local compressive stress, 
which may provide a driving force for nanoline nucle- 
ation. We demonstrated earlier (Fig[7J) that the 7-5-7 
structural motif may be achieved by a simple bond rota- 
tion process. Using this motif, the proposed mechanism 
suggests that the central four atoms between the two Bi 
dimers in the proposed mechanism twist, so that a zigzag 
core structure is formed. From this the Haiku structure 
evolves, with a maximum barrier in the calculation of 
less than 1.5 eV. The reaction pathway and associated 
energy diagram are shown in Fig. 1201 This mechanism 
has a number of interesting features, but also a major 
limitation. The calculation takes place in a periodically 
repeated unit cell one dimer row wide, and is therefore 
considering the nuclcation of an infinite nanoline, with- 
out ends. It is suggested that a long nucleus forms by the 
filling of several unit cells of a missing dimer trench, but 
this is an unstable situation |129| . Moreover, in our STM 
observations, we observe the appearance of very short 
nanolines, which then grow. The termination of these 
finite nuclei must therefore be considered. STM obser- 
vations of nanoline evaporation show that nanolines do 
not break up in the middle, instead they evaporate from 
the ends. From this, it may be inferred that the energy 
of termination is quite high. While there has been no 
systematic study of termination, we have performed cal- 
culations on many and varied structures, all of which 
have been found to have a high energy. In tightbinding 
modelling, the energies were at least 1 eV/end , while 
those terminations which we have considered for the pro- 
posed zigzag nucleus structure have much higher ener- 
gies, around -1-1.6 eV/end, giving a total nucleus energy 
of over +A eV. This suggests both that understanding 
the termination of the nanolines is an important area for 
future research, and that nucleation structures should be 
considered in 2D, rather than in ID (which may not be 
feasible with a technique more accurate than tightbind- 
ing). Given that the prevalent structures on the Bi-rich 
surface are Bi surface dimer/ IDV complexes, alternative 
nucleation mechanisms which use these structures as a 
basis may be more fruitful and are being explored. 

4-.3 Electronic Structure 

The knowledge of the electronic structure of the Bi nano- 
line has stemmed from two sources: variable- voltage STMj2j 
I132lll3,'i| and DFT calculations of the electronic structure |^ 

in terms of the local density of states (LDOS), or the 
Tersoff-Hamann approach to STM, which uses the LDOS). 
Knowledge of the atomic structure of the Bi nanoline al- 
lows detailed calculations of its electronic structure to 
be made, which can then be compared to the observa- 
tions made in STM. While the majority of the electronic 
contrast is expected to be dominated by the Bi dimers 



Fig. 21 10 nm x 10 nm images of a Bi nanoline. The sample 
bias voltages used are -2.0 V, -1.2 V, -0.6 V, -0.3 V, -1-1.5 V 
and -0.8 V, in (a-f) respectively. As the sample bias is re- 
duced, between (a) and (d), the nanoline changes contrast 
from light to dark relative to the surrounding Si(OOl). Over 
this range, some of the dimers in the nanoline (marked A in 
(c) and (e)) exhibit a different voltage contrast. At very low 
biases, around -0.3 V, an enhancement of the dimers around 
the nanoline, similar to that seen around a missing dimer de- 
fect in clean Si(OOl) '1481, is seen. This is visible (and marked 
schematically as B) in (d). In (f), the resolution is sufficient 
to see that the corrugation of the Si dimers closest to the 
nanoline is increased, (marked schematically by the dotted 
black lines) suggesting a greater separation, and hence ten- 
sile strain. Reprinted from Surf. Sci. Lett. 527, Owen et al. 
Interaction between strain and electronic structure in the Bi 
nanoline, p. 177, (2003), with permission from Elsevier. 

in the nanoline, other details of the LDOS have allowed 
a stronger identification of the Haiku structure with the 
Bi nanoline to be made, with the other possible candi- 
date structures ruled out after comparison with the STM 
data. 

We begin by summarizing the data available from 
STM, which is illustrated in Fig. EH 

1. At high biases, the Bi nanolines appear bright com- 
pared to the surrounding Si(OOl) surface , 1 36. 1 46iJ^ 

2. At low biases, the Bi nanolines become dark com- 
pared to the surrounding Si(OOl) surface 9- 125.. 1251 

3. At low biases, the Si dimers neighbouring the Bi 
nanoline show enhancement, or brightening |132| sim- 
ilar to that seen around missing dimer defects on 

Si(ooi)[TlHIIISII. 

These facts suggest various things: first, that the 
Bi dimers are physically higher than the surrounding 
Si(OOl) surface dimers, as at high biases the contrast 
in STM is determined primarily by geometrical struc- 
ture; second, that the local electronic structure of the 
Bi dimers will show states further from the Fermi level 
than the surrounding Si dimers 1 32? 1 3311 1 38| . as this de- 
termines the STM current at low biases; finally, that the 
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Bi nanoline structure must be under tensile stress, which 
causes the Si diniers near the nanoline to be pulled to- 
gether, away from their relaxed surface structure, which 
is the origin of enhancement at low biases [TlSlllSOj . 

4-.3.1 Electronic Properties of the Bi dimer Examina- 
tion of the detailed electronic structure of the Miki and 
Naitoh models |T!?D, 131 , 55 and the Haiku structure ji;-{2l 
I14(JII133| all show that the Bi dimers become dark at 
low bias voltages by contrast to the surrounding sili- 
con dimers. The calculated DOS for the Miki and Haiku 
structures are shown in Figure [521 138 • The vl peak 
comes from the up atoms of the substrate Si dimers, 
as does the v2, while the v3 peak comes from both the 
Bi and Si dimer bonds and vA from the Bi atoms in the 
line models. When the substrate dimers are passivated 
with H, the peaks vl and v2 almost disappear leading 
to peaks in the gap which arise from the Bi dimers (v3 
and vA). This data confirms earlier modelling and bias- 
dependent STM measurements [TS^ discussed below. 

Consideration of the electronic structure of the Bi 
dimer reveals that it is rather unusual. In its native 
form, Bi forms puckered hexagonal sheets, similar to the 
structure taken by As. Full-potential LAPW modelling 
of Bi in a bulk diamond cell (thus forcing it to take 
up tetrahedral bonding) showed that there is very little 
hybridisation between the s- and D-orbitals |128j . Thus 
we expect Bi to tend towards p'^ bonding with a filled 
s-orbital (as a lone pair) and to be energetically most 
favourable with 90° bond angles. This explains to some 
extent the Haiku structure, where bond angles are al- 
most 90°. Projection of the local densities of states calcu- 
lated using pseudopotential-based DFT calculations 
suggests that the Bi dimer in the Haiku structure has hy- 
bridisation between s and Pz, with px and Py separated 
off. The FLAPW calculations showed that the relativis- 
tic effects of the full core increased the s-p splitting in 
the Bi atom; while relativistic effects can be incorporated 
in pseudopotentials, the Haiku structure might be more 
stable, and show bonding closer to p^, if these effects 
were taken fully into account^. 

The Bi-Bi bond is a rather weak bond, with the Bi-Si 
bonds providing most of the energetic stability. This can 
be seen by considering the electronic structure of a Bi 
ad-dimer on the clean surface^. The electron localisation 
function (ELFl fTM] for this ad-dimer is shown in Fig- 
ure l23r al |152| . The ELF gives a quantitative analysis of 
the bonding properties of a system, with a value of 0.5 
equivalent to the localisation seen in the homogeneous 
electron gas and a value of 1.0 equivalent to complete lo- 
calisation. The figure shows an isosurface with ELF=0.8 

We note that a full potential calculation on a unit cell 
large enough to hold the haiku structure is probably compu- 
tationally prohibitive 

While not a perfect mimic of the Bi dimer in the Haiku 
structure, the smaller cell allows more detailed calculations 
to be made. 
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Fig. 22 The calculated DOS of the Miki and Haiki struc- 
tures on a clean surface (top) and a H-terminated surface 
(bottom) are shown. The marked states vl to v4 are de- 
scribed in the text. Reprinted with permission from Ref 
Copyright (2005) TOP Publishing. 



(equivalent to reasonably strong covalent bonds). It can 
be seen that the Si-Si bonds are rather strong, as ex- 
pected, and that the Bi-Si bonds are still visible at this 
level, though they are weaker than the Si-Si bonds. The 
Bi-Bi bond is not visible (in fact it has a maximum value 
around 0.7) but the Bi lone pairs are visible as rather dif- 
fuse, spherical objects (by comparison to the up atoms of 
the substrate dimers). This implies that the Bi-Bi bond 
is rather weak, and the lone pairs would not be avail- 
able for bonding 152' . The energy barrier to breaking 
the bond in this Bi dimer has been calculated with DFT 
to be 0.17eV going from over the dimer row to over the 
trench between rows and 0.15 eV going in the other di- 
rection, and the barrier is shown in Fig. I23r b). Having 
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the Bi dimer over the dimer row is more stable than over 
the trench between rows by about 0.02 eV. This imphes 
that there should be a finite chance of imaging the Haiku 
structure with the Bi dimers out-of-phase with the un- 
derlying Si dimer rows. However, we note that with the 
energy difference noted above, the dimer would be over 
the row 98% of the time at room temperature, and the 
likelihood that a whole line of dimers would flip, will 
be even smaller. There are occasional images with the 
dimers out-of-phase, but these are rare. 




Fig. 23 (a) The electron localisation function (ELF) of a 
Bi dimer on Si(001)[ni| for a value of 0.8. The ELF indi- 
cates the localisation of pairs of electrons, indicating strong 
bonding and lone pairs. The figure shows the weakness of 
the Bi-Bi bond, and the diffuseness of the Bi lone pair. See 
text for full discussion, (b) Energy barrier for Bi dimer to go 
from over dimer row to over trench between dimer rows. The 
dotted line is a spline fit intended as a guide for the eye. 

In the next section, details of some elements which 
will attack the Bi-Bi bond will be presented, along with 
others that will not. Here we consider the local structure 
of the bond. When the Bi-Bi bond is over the Haiku core, 
it is almost perfectly relaxed, with bond angles extremely 
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Fig. 24 A comparison of simulated STM for the 3 proposed 
models. An STM image of the Bi nanoline has been aligned 
with the centre of the calculated nanoline for each structure. 
For the Miki model, the background Si dimers are offset rel- 
ative to the experimental image, while for the Naitoh model, 
the Bi dimer spacing is much too small. For the Haiku model, 
both of these aspects fit very well. The change in contrast of 
the Bi nanoline, and the dimers adjacent to it, at low biases 
is discussed in the text. The data for the Miki and Naitoh 
models are reprinted with permission from Ref |55) . Copy- 
right (2004) by the American Physical Society. 

close to 90°. Insertion of any species into the Bi-Bi bond 
will distort these angles (and the angles formed with 
the underlying Si atoms), with a potentially large en- 
ergy penalty. Thus the Bi-X-Bi angle should be relatively 
close to 180° with short Bi-X bond lengths if an element 
X is to insert successfully into the bond. Furthermore, 
the lone pair on the Bi dimer, which might be expected 
to be a strong electron donor, is in fact rather passive, 
as it lies some distance below the Fermi level [152] . The 
implications of this for the chemical reactivity of the Bi 
nanoline are discussed below in Section ^31 

4-3.2 Strain-induced Enhancement The phenomenon of 
strain-enduced enhancement at low STM bias voltages 
has been explored before [n5lll50| : we reproduce the ar- 
guments here for convenience |132| . The Si(OOl) surface 
shows a dimer reconstruction which provides three bonds 
per surface Si atom (compared to two bonds per atom 
for the bulk terminated surface), but which pulls the 
bond angles and distances away from the ideal tetrahe- 
dral values. The net result is that the local band gap 
is reduced (in simple chemical terms, the hybridisation 
towards sp"^ is decreased, leading to a decreased split- 
ting between the resultant hybrid bonds). If the surface 
dimers are then further distorted, their local states will 
move closer to the Fermi level. The single missing dimer 
defect (IDV) has second layer atoms bonding across the 
defect which distorts the dimers either side of the defect 
(though lowering the energy because the second layer 
atoms have no dangling bonds). This distortion in turn 
causes the states near the Fermi level on the neighbour- 
ing dimers to move towards the Fermi level, leading to 
a brightening or enhancement at low biases. 

Low-bias simulated images of the Naitoh and Miki 
structures | 131b130>i55| show that the Si dimers adjacent 
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to the Naitoh become darker than the surrounding Si 
dimers as the Fermi level is approached, while those 
on the Miki model are unchanged. By contrast, with 
the Haiku structure, the Si dimers immediately adja- 
cent to the nanolinc become brighter closer to the Fermi 
level |l.S2t..l40 ■ in good agreement with the STM images. 
A comparison of high and low-bias simulated images of 
the 3 structures are shown in Fig. |^ However, the pre- 
diction from the Miki model does agree very well with 
the appearance suggested from STM images [51ll49j . in 
which the shoulder dimers around a IDV on a Bi-rich 
surface become dark at low bias voltages. On the basis 
of these calculations, we suggest that the lack of enhance- 
ment at low bias voltages indicates almost complete re- 
laxation of the strain of the IDV by the presence of the 
Bi dimers in the Miki model. 

Taking the principle of strain-induced enhancement 
further, the local subsurface strain could be determined 
by cross-sectional STM of the Bi nanoline, of which there 
is none, or from LDOS calculations, by projection of the 
charge density on a plane perpendicular to the (001) sur- 
face. We show the projected charge density from GGA 
calculations for the clean Si(OOl) surface, the Haiku struc- 
ture and the Haiku structure with hydrogen on the Si(OOl) 
surface in Fig. ESl 132, . on the (110) plane, parallel to a Si 
dimer row. By comparing in particular the clean Si(OOl) 
surface and the clean Haiku structure, the change in dif- 
ferent states induced by the Haiku can be seen. Begin- 
ning with the clean surface, the localisation of the states 
close to the Fermi level in the top few surface layers 
shows the strain induced by dimerisation. In the Haiku 
structure, these states are even more strongly localised 
on the dimers immediately adjacent to the Haiku struc- 
ture, which causes their relative enhancement in STM. 
Meanwhile the triangular Si substructure of the Haiku, 
c.f Fig. |S1 is somewhat relaxed compared to the clean 
Si dimers, although it is strained relative to the hydro- 
genated dimers and perfect bulk Si. This goes some way 
to explaining the stability of the subsurface 5-7-5 ring 
structures, which serve as a highly effective relief mech- 
anism for the epitaxial stress exerted by adsorbed Bi. 
The absence of states close to the Fermi level suggests 
that the nanolines are likely to block surface conduction 
perpendicular to the nanoline, and are not likely to act 
as a nanowire. 

The lack of any states on the Bi near the Fermi level 
is clearly the cause of the darkening of the nanoline rel- 
ative to the Si(OOl) surface seen in STM at low volt- 
ages. In the full charge densities (Fig. I^SJiv)), the den- 
sity of states associated with the Bi dimers, combined 
with their higher physical height, explains their relative 
brightness in STM at higher bias voltages. In the case of 
the hydrogenated surface, the Bi dimers have a similar 
charge density as in the clean surface, but in this case, 
the Si vr-bonds have been eliminated, so the Bi dimer 
remains bright in STM at all biases [T^ . 



4-. 4 Chemical Properties 

As described above in more detail, the bonding orbitals 
in the Bi dimers have essentially a p'^ bonding charac- 
ter, while the lone pair is a diffuse s-type orbital, some 
10 eV below the Fermi level 128 . Thus the Bi dimer is in 
a very stable electronic state, and is not expected to be 
very reactive. Moreover, this low-energy configuration is 
dependent upon the maintenance of bond angles close to 
90° , which will have an effect on the energetics of any in- 
sertion into the Bi-Bi dimer bond. However, the low-bias 
images shown above reveal that although the nanoline 
appears clean at high biases, some contrast between dif- 
ferent dimers is visible at lower biases, suggesting that 
some chemical attack has occurred. 

4-4- 1 Reactivity of the Bi dimer We have proposed that 
the Bi nanolines, although semiconducting themselves, 
might be used as templates for deposition of other ma- 
terials, such as metals, active molecules, nanoparticles 
etc. |135lll53| . In order for the nanoline to have utility 
in this role, it must be possible to mask the substrate 
around the nanoline, so that adsorbants stick preferen- 
tially to the nanoline. It is also important, for any nano- 
electronics applications, to have some method for ren- 
dering the substrate insulating after growth of the nano- 
line. The chemical reactivity of the Bi dimers, and the 
nanoline as a whole, is therefore crucial to any implemen- 
tation of this concept. Atomic hydrogen has been used 
as a mask on the bare Si(OOl) surface for STM lithogra- 
phy exDeriments pil79ll80| . in which an STM tip is used 
to remove lines or areas of hydrogen from the surface, 
and thus provide areas onto which metal will adsorb, 
as in Fig. El However, the writing process is slow, and 
lacks scalability. Use of the Bi nanolines as templates 
would remove the need for this writing, while providing 
long, and more uniform patterns. Hydrogen is therefore 
a natural species to use as a mask around a Bi nanoline. 

One of the most favourable properties of silicon as 
an electronics material is the high quality of its oxide 
film. Oxidation of the silicon around the Bi nanoline is 
therefore a natural way in which to isolate the nanoline 
electrically from the substrate. The usual method of ox- 
idation of silicon is by exposure to molecular oxygen or 
to water at elevated temperatures, but this only gives 
microelectronics-quality oxide when grown at ca.8GG°C. 
Since the Bi nanoline is not stable above 600° C, ozone 
was considered as an alternative for low-temperature ox- 
idation. Ozone preferentially attacks the backbonds of 
the surface Si dimers, forming a stable oxide, even at 
room temperature'154 . More recently, the need for high- 
k dielectrics has driven the development of silicon nitride 
and oxynitride as alternatives to silicon dioxide. Nitrida- 
tion is usually done using ammonia, which breaks up on 
contact with the clean silicon surface above 200 K |155| . 
As an added bonus, when ammonia reacts with silicon, 
the nitrogen moves below the surface, but the hydrogen 
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Fig. 25 Contour plots of charge density for the states (i) 0.2 eV, (ii) 0.6 eV and (iii) 1.0 eV below Ef, and (iv) the total 
charge density for: (a) the clean Si(OOl) surface (first column); (b) the Haiku structure on the clean surface (second column), 
and (c) the hydrogenated surface (third column). The 0.2 eV image for the hydrogenated surface has been replaced with a 
ball-and-stick model of the Haiku structure, as there are no states within 0.2 eV of E_f. (The highest electron density contours 
are 0.003, 0.060, 0.070 and 0.500 electrons/cubic Ain the lst-4th row respectively of columns a,b. The maximum contours in 
column c are 0.01,0.05,0.50 electrons/cubic Ain the 2nd-4th row.) See text for a full discussion. Reprinted from Surf. Sci. Lett. 
527, Owen et al. Interaction between strain and electronic structure in the Bi nanoline, p. 177, (2003), with permission from 
Elsevier. 



remains on the surface, passivating it. It can therefore 
play a dual role, as insulating reagent, and as masking 
material. The reaction of these species with the Bi nano- 
line are discussed in the sections below. 

4-4-^ Hydrogen The Bi nanoline has been shown to be 
almost completely inert to attack by atomic hvdrogen |126l 
ISMUni- FigEHa) shows a close -up image of a nanoline 
on a Si(001):H monohydride surface. In this image, nano- 
lines were formed at 590°C, and the surface was cooled 
to 330°C, the surface was saturated with atomic H. The 
surface has the monohydride phase, with the Si-Si a- 
bond intact, and the H atoms saturating the dangling 
bonds. Thus each Si-H dimer is resolved as a pair of 
dots. Some white dots on the Si surface may be holes 
in the monohydride layer, stray bismuth atoms, or ad- 
sorbed -OH resulting from water impurity in the hydro- 
gen. While the nanoline is almost completely unreacted, 
one dimer has become dark. This sort of defect is ex- 



tremely rare. The preferential adsorption of the H onto 
the Si is consistent with DFT-LDA modelling, which 
found that H on a Bi dimer is 3 eV worse than H on 
a Si dimer. There is therefore a strong driving force for 
any H which does adsorb onto the Bi nanoline to diffuse 
off it, and onto the Si. 

4-4-3 Oxygen The oxidation properties of a monolayer 
of Bi on Si(OOl) shows an unusual changeover in be- 
haviour with oxygen exposure. An EELS/ Auger study 
found that for small exposures, the reactivity of a Bi 
overlayer on Si(OOl) is 100 times less than the reactiv- 
ity of oxygen with the bare silicon surface. However, 
for large exposures, above ca. 3 x 10^ L O2, the oxygen 
Auger intensity continues to increase for a Bi/Si(001) 
surface, at a point where the signal for the clean Si(OOl) 
surface has saturated. Thus for very large exposures, the 
Bi promotes oxidation of the silicon. Moreover, the stoi- 
chiometry of the oxide formed on the Bi/Si(001) surface 
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Fig. 26 a: A fiUed-states atomic-resolution image of a H- 
terminated surface, showing a section of a nanoline. The 
imaging conditions here are -2.5 V, 0.1 nA. (b); Bi nano- 
lines after exposure to 2L O3 at 300 K. There is virtually 
no attack on the nanolines, while the substrate has been oxi- 
dised, (c): After exposure to 150L O2 at 700 K and an anneal 
to 840 K, some gaps appear in the nanoline, but most of it 
is unaffected, (d): Exposure of a H-terminated surface with 
Bi nanolines to air results in complete oxidation of the nano- 
lines, which now show up dark relative to the H-terminated 
Si. 



nealing this surface introduced some defects into the 
nanoline. After annealing at 450°C, gaps appear in the 
nanoline, which comprise 1-2 unit cells. This trend con- 
tinued with a further anneal to 570°C, as shown in Fig. 
I26r c'). Such piecemeal damage may be explained by the 
observation that the activation barrier for Bi to desorb 
from oxidised Si has been found to be 0.7 eV lower than 
for unoxidised Si |156j . The small contrast changes, which 
were seen after ozone exposure at room temperature 
may therefore indicate reaction with Bi-Si, or Si-Si back- 
bonds around the nanoline, so making that unit cell more 
susceptible to evaporation during high-temperature an- 
neaUng. Molecular oxygen, when adsorbed at 400°C, be- 
haved in a very similar way to ozone. Again there was 
preferential attack of the silicon, leaving the nanoline 
mostly intact, but also gaps appeared in the Bi nanoline 
after annealing for about one hour'135'. 

While the Bi nanoline is remarkably resistant to ox- 
idation in UHV, it quickly oxidises when exposed to 
air, as shown by the following simple experiment. A H- 
terminated Si(OOl) surface with Bi nanolines was taken 
out of the imaging chamber into the UHV system load- 
lock, which was then vented to atmosphere for 30s, and 
then immediately pumped back down to UHV. The re- 
sulting surface is shown in Fig. I26f d). The nanolines, 
which previously were bright at all biases relative to the 
H-tcrminated background, are now black lines, indicat- 
ing that they have suffered significant attack. This agrees 
well with the previous study, which found that for very 
large exposures to oxygen, the Bi dimcrs promoted oxi- 
dation of the underlying silicon. 



is closer to Si02, which would imply better insulating 
properties. Semiempirical calculations found that inser- 
tion of oxygen into the Bi-Bi bond is energetically much 
less favourable than for Sb or As (0.32 eV vs ca. 4.7 eV). 
Moreover, insertion into a Bi-Si bond or a Si-Si bond is 
much more favourable (4.5 eV and 5.5 eV respectively), 
so that preferential attack of the silicon is expected. 

The reactions of molecular oxygen and ozone with 
the Bi nanoline show a similar changeover in behaviour |135j 
The surface was exposed to ozone at room temperature 
and, for comparison, to molecular oxygen at 400°C. Both 
surfaces were annealed up to 570°C, near the stability 
limit of the Bi nanoline. The room-temperature ozone 
experiments, using exposures up to 20 L of O3, revealed 
that the silicon background could be saturated without 
visible damage to the nanoline. In large-scale images, 
the Bi nanoline appears to be completely unaffected, as 
shown in Fig.jSBfb). However, in small-scale images, and 
by varying the imaging bias, some subtle differences be- 
tween the appearance of various dimers in the nanoline 
could be seen |135j . This variation may be due to attack 
of Bi-Si backbonds, as found to be favourable in MNDO 
simulations^^, which would leave the Bi-Bi dimer in- 
tact, but perturb its electronic structure somewhat. An- 



4-4-4 Ammonia The adsorption of ammonia onto Si(OOl) 
as a precursor to silicon nitride or oxynitride growth has 
been studied by a variety of experimental and theoretical 
techniques, which concluded that at room temperature, 
the ammonia dissociated on the surface, forming NH2 
and H groups |TK5| . We have found that ammonia at- 
tacks the silicon preferentially at room temperature. An 
image of the Bi:Si(001) surface with ammonia termina- 
tion is shown in Fig l27l Within each dimer, the light grey 
end has been identifed with the NH2 group, while the 
dark grey end is the H atom |15S] . The apparent height 
difference in STM matches the physical height difference 
between these two species. There is some ordering of the 
NH2 groups, such that zigzag patterns and straight lines 
are seen running along the dimer rows. 

Reaction of ammonia with the Bi nanoline might be 
expected to proceed by the initial formation of a hydro- 
gen bond between an ammonia H atom and the Bi lone 
pair, followed by insertion into a Bi-Bi or Bi-Si bond. 
However, as discussed above, the lone pair on the Bi 
dimer is essentially inert [152] . During experiments on a 
surface which was pre-exposed to atomic H to saturate 
the Si dangling bonds, for exposures of up to 70L of am- 
monia at substrate temperatures of up to 500K, the Bi 
nanoline was largely unaffected. However, adsorption of 
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Fig. 27 A 64 nmx64 nm STM image of tlie ammonia- 
terminated Bi;Si(001) surface. Tiie background Si is com- 
pletely terminated by a mixture of H and NH2 groups. Dark 
patches on the Bi nanoline reveal some attack by ammonia. 
The inset shows that reaction occures within the space of a 
single dimer, and the resulting feature is asymmetric, much 
like the termination of the background Si. 

ammonia onto a Bi nanoline surface where the Si dimers 
are not passivated does result in some reaction with the 
nanoline, even at room temperature. This difference may 
indicate that ammonia reacts first with the silicon, form- 
ing active NH2 species, which are mobile and can attack 
the Bi dimers. Calculations suggest that it is energeti- 
cally favourable for an NH^, group to insert into the Bi-Bi 
dimer bond, and form a symmetrical Bi-NH-Bi feature. 
However, STM images of reacted Bi dimers shows an 
asymmetric feature, not unlike the feature seen on back- 
ground Si dimers. Examples may be seen in the inset to 
Fig. Nevertheless, the probability of attack by am- 
monia is low, providing the promise that ammonia can 
be used to grow silicon nitride without damage to the 
nanoline. 

4.4-5 Metals Metal deposition onto the Bi nanolines 
is currently under investigation. Deposition of Ag onto 
Bi nanolines on the clean Si(OOl) surface 159' reveals 
a strong preference for adsorption onto the Si, leaving 
the Bi nanoline clean (in complete contrast to the be- 
haviour of ErSi2 NWs on Si(OOl) discussed above in 
Sec. I3.3|l . However, the passivity of the Bi nanolines to 
atomic hydrogen and ammonia allows these species to 
be used as nanoscale masks. Ag is known not to ad- 
sorb onto H-terminated Si(001)|Hn| as shown in Fig. |21 
while In completely dewets from the Si(OOl) surface, 
forming large droplets when H is adsorbed onto an In- 



covered surface pi^il) . In all cases, following metal de- 
position onto an ammonia-terminated Si(OOl) with Bi 
nanolines, preferential adsorption onto the Bi nanoline 
(or onto background Bi) was seen, indicating that the 
ammonia-terminated surface is very effective as a mask 
against metal deposition. The reaction of the metal with 
the nanoline template is more complex, however. Two 
types of behaviour have been seen with different met- 
als, which may be described as wetting and non-wetting 
behaviour. Wetting occurs when the interaction of the 
metal with the nanoline is stronger than the interaction 
of the metal with itself. Nonwetting occurs when the in- 
teraction of the metal with itself is stronger. Examples 
of each type of behaviour for different metals. In and 
Ag, are shown in Fig. [2HI Adsorption of In results in 
significant intermixing of In and Bi, forming a zigzag is- 
land structure, which is thought to result from a chain 
of In and Bi atoms, as shown in the inset to Fig. I28r a). 
Second-layer islands have a distinctive bright hexagonal 
feature, which is 19 A or 2.5 dimer rows wide. This ra- 
tional relationship between the In island dimension and 
the underlying Si unit cell distance demonstrates that 
these In islands grow epitaxially on the nanoline tem- 
plate, raising the possibility of the growth of long single- 
crystal nanowircs on the Si(OOl) surface. By contrast, 
adsorption of Ag results in the formation of small clus- 
ters of Ag, around 0.5 nm in height, at a very early stage 
of deposition, as shown in Fig.|2HIc). Further deposition 
at room temperature results in an increase of the num- 
ber density of these nanoclusters, but the modal size 
increases only slightly from 0.5 nm to 0.6 nm, with a 
strongly peaked height distribution. This behaviour may 
be thought of as non-wetting behaviour. 

4-5 Burial 

While the Bi nanoline has a number of interesting prop- 
erties, for any nano-electronics device application it must 
be possible to passivate them against atmospheric at- 
tack. One way in which to do this is to bury them in 
a layer of epitaxial silicon. Moreover, the burial of the 
Bi atoms may lead to their acting as a 2D delta-doping 
layer, with close to atomic-layer precision in the depth 
placement. However, as with other nanostructurcs, for 
example burial of InAs/GaAs quantum dots in GaAs, 
the process of burial is likely to change the shape and 
detailed structure 

Simple overgrowth of a surface containing Bi nano- 
lines results in the destruction of the nanolines ' 1 6 1',' 1 42] . 
As shown in the Secondary Ion Mass Spectrometry(SIMS) 
profile in Fig. I29r a). after overgrowth of Si the vast ma- 
jority of the Bi has segregated to the surface of the grown 
film, although there is a small peak in the Bi concentra- 
tion at the original surface position. Furthermore, X- 
ray analysis'1421 shows that the 1-D character of the 
nanoline is lost after capping, even by an amorphous 
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Fig. 28 (a): Indium reacts with the Bi nanoline to form a 
zigzag feature, (b): A two-layer In island shows hexagonal 
features, which are 19A or 2.5 dimer rows wide, revealing 
an epitaxial relationship with the Bi nanoline. (c): Isolated 
0.5 nm Ag nanoclusters form after deposition of Ag atoms at 
room temperature onto the Bi nanoline surface, (d): Height 
Analysis reveals a strongly-peaked size distribution, with the 
peak around 0.5 nm, which may correspond to a magic cluster 
of 13 atoms. 



overlayer deposited at room temperature. The mecha- 
nism by which the Bi floats to the surface during growth 
is thought to involve exchange of subsurface Bi atoms 
with surface Si atoms, as Bi is an effective surfactant 
during growth. In order to prevent this, a further sur- 
factant layer of Bi can be deposited around and over 
the Bi nanolines, before Si overgrowth' 161) . In this case, 
the exchange mechanism is inoperative, as the second- 
layer Bi in the nanolines will not exchange with first- 
layer Bi, and so it is possible to block the segregation 
of the Bi nanolines to the growth surface. In the SIMS 
profile in Fig. I29r b). a thin crystalline Si overlayer has 
been grown at 400°C after deposition of a Bi layer, and 
this has been capped by depositing further Si at room 
temperature, so as to bury the growth surface in amor- 
phous Si. In this case there, are two peaks in the SIMS 
profile. 'A' refers to the Bi density at the original Bi 
nanoline surface, while 'B' refers to the surface of the 
crystalline overgrowth. The surface of the overlayer has 
approximately 1 ML of Bi on it, from which the quantity 
of buried Bi can be determined. This is approximately 
5% ML, which is typical for a Bi nanoline surface. Thus, 
using the Bi surfactant layer, most of the Bi nanolines 
have been preserved at the original surface. In order to 
gain more information about the exchange process, a 
surfactant layer of Sb was also used, instead of the sur- 
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Fig. 29 (a): A SIMS profile of a Bi nanoline surface, formed 
at the position marked 'A' and then overgrown with Si at a 
substrate temperature of 400° C. Most of the Bi has segre- 
gated to the surface, and the nanolines have been destroyed, 
(b): The Bi nanoline surface (A) has been capped with a 
monolayer of Bi at 400°C, before overgrowth. After over- 
growth at 400°C up to position 'B', the growth layer was 
buried in amorphous Si by further deposition of Si at room 
temperature. In this case, nearly all the Bi from the Bi nano- 
lines remains at the original surface position, 'A', (c): As (b), 
except that the capping layer material was Sb, in place of Bi. 
Here most of the Bi has segregated to the growth surface (B), 
leaving Sb at the original surface (C), indicating that the Bi 
has exchanged with Sb from the capping layer. 
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Fig. 30 (a): A ball-and-stick model of a proposed struc- 
ture for a buried Bi nanoline. The haiku structure has recon- 
structed into a 2DV-hke trench, with the Bi dimers attached 
to Si overhead, (b) : A cross-sectional TEM image of Bi nano- 
lines buried in a Si epitaxial layer. The one on the left is defect 
free, while the one on the right has triggered dislocations. 



factant layer of Bi. The resulting SIMS profile is shown 
in Fig. I29r c). In this case, there is a large peak of Sb at 
the original nanoline surface position 'B', while a con- 
siderable fraction of the Bi has segregated to the surface 
of the overlayer 'C. This result indicates that there is 
considerable exchange between the Sb surfactant layer 
and the Bi nanolines, and by this method, a layer of Sb 
has been successfully buried, with a sharp peak at the 
original nanoline surface. 

While this exchange process is successful at burying 
the Bi atoms, it does not provide any information about 
the structure of the buried Bi. For device applications, 
it is important to retain the 1-D character of the orig- 
inal nanoline. It is also of interest whether the Haiku 
structure will be preserved. STM cannot observe sub- 
surface features, and instead the structure of the buried 
Bi nanolines was studied using an X-ray method |l62) . 
The results of this analvsis |l41| indicate that two im- 
portant elements of the Bi nanoline structure have been 
preserved: the 1-D character, and the dimerisation of 
the Bi atoms. Tightbinding and DFT calculations of 
the buried Bi nanoline structure found that the Haiku 
structure was not stable, and broke down into a hol- 
low 1-D structure with Bi dimers bonding to the sili- 
con grown over the nanoline, and the Si substructure 
reconstructing into something similar to a 2DV surface 



defect. A ball-and-stick model of this proposed struc- 
ture is shown in Fig. I^Uf a). It has significantly lower 
energy per Bi atom than other proposed structures, in- 
cluding substitutional and interstitial Bi defects. This 
structure may be compared to the cross-sectional TEM 
image of buried Bi nanolines in Fig. l^Of b). Two nanolines 
are present within the TEM image. The one on the left 
has strain associated with it, and a dark core, which may 
match with the proposed structure. The one on the right 
has triggered misfit dislocations along {111} planes, c.f. 
Fig. IHl These dislocations, which indicate the presence 
of a large local stress, may suggest that the Haiku sub- 
structure of the nanoline has survived overgrowth. A re- 
cent X-ray Standing Wave (XSW) study of a Bi nano- 
line sample which had been capped with amorphous Si 
at room tempcrature |144| has cast doubt on the validity 
of the Haiku model. However, more detailed analysis of 
lines buried with amorphous silicon n^T| showed that the 
dimerisation and one dimensional character of the nano- 
lines is destroyed by burial in amorphous silicon (and, 
as discussed above, burial in crystalline silicon changes 
the Haiku reconstruction while leaving the one dimen- 
sional, dimerised character of the nanolines unaffected). 
Furthermore, the presence of a significant quantity of Bi 
at the growth surface may have an effect on the results 
obtained, unless it is removed |142 |. 

5 Conclusions 

We have given a broad overview of the different mate- 
rial systems and experimental methods which result in 
self-assembled ID nanostructures on semiconductor sur- 
faces, particularly the technologically important Si(OOl) 
surface. There are many different self-organising systems 
which result in ID nanoscale features on the surface of 
Si(OOl) and Si(lll). Aside from the Pt/Ge(001) system, 
none of these approach the length, and degree of per- 
fection, of the rare-earth silicide wires and the Bi nano- 
lines. The metal/Si(lll) nanowire systems will tend to 
form 3 equivalent domains, limiting the long-range or- 
der of the nanowires, while the methods which rely on 
step-edge adsorption are limited by the spacing of step 
kinks (though this can be controlled with careful sam- 
ple preparation). The rare-earth silicide nanowire family 
are not nanowires per se (in that they are a metastable 
state, not seen in a bulk structure), but are better re- 
garded as a conventional heteroepitaxial system in which 
the lattice mismatch is anisotropic. Even so, while many 
different rare earths superficially appear to behave in a 
similar fashion, review of the literature indicates that in 
fact several different nucleation and growth mechanisms 
are responsible for the formation of nanowires within 
this family, both kinetic and thermodynamic. The recent 
observation of nanowires with an orthorhombic crystal 
structure 11 7J, rather than the hexagonal structure pre- 
viously reported, suggests that there is much more to be 



Self-assembled nanowires on semiconductor surfaces 



29 



learnt about this family of materials. Experimental stud- 
ies of all these systems has so far been largely confined 
to STM, apart from X-ray and TEM studies of buried 
Bi nanolines, and TEM of the endotaxial RE nanowires. 
A much greater understanding of the physical, chemical 
and electronic properties of these systems, would come 
from other experimental techniques, such as optical and 
vibrational spectroscopic techniques. (For example, the 
strained Haiku core of the Bi nanoline may have a sig- 
nature bond vibration frequency.) 

The main focus of this review has been the Bi nano- 
line system. We have described the physical and elec- 
tronic structure of the Bi nanoline, its reactivity with 
different materials, and the surface phenomena associ- 
ated with it. Although there has been no direct observa- 
tion of the core structure of the Haiku model, it provides 
a natural explanation for all the observed experimental 
phenomena, and is well-supported by experimental and 
theoretical data. The fundamental source of many of the 
properties of the Bi nanoline is the subsurface core of Si 
which has recrystallized into a hexagonal phase, Lons- 
daleite (known from meteorites). The formation mech- 
anism of this core structure remains unknown. A sig- 
nificant activation barrier is expected, given the high 
temperature and long incubation time required to form 
the nanolines, and the co-existence of the nanolines with 
the Miki structure, which is less stable, but kinetically 
easy to form. While mechanisms proposed thus far have 
concentrated on the twisting of pairs of atoms, an alter- 
native may well be a stacking- fault- like shift of the entire 
core from the cubic phase to the hexagonal phase. Given 
their structure, it is clear that the Bi nanolines are not a 
conventional heteroepitaxial system like the silicide fam- 
ily. There is no variation in width, either with annealing 
temperature or time - only the number density of nano- 
lines and the length changes with further deposition and 
annealing. The nanolines are only ever one unit cell wide, 
although clusters of nanolines are sometimes seen. The 
nanolines are not the result of simple adsorption onto 
symmetry-breaking features of the surface, such as step 
edges, or linear defects. Nor are they a periodic recon- 
struction of the surface, as in the Pt/Ge wiresj^]. They 
might best be described as an adsorbate-stabilized sur- 
face recrystallization. 

The chemical properties of the nanoline stem from 
the behaviour of the Bi dimer with its preference for 
p^ bonding and an inert s-type lone pair, rather than 
the sp'^ hybridization of Si dimers. Although the Bi-Bi 
bond is weak, it is also passive, and only In shows any 
significant reaction with it, probably because the zigzag 
structure produced allows the Bi atoms to retain their 
preferred 90° bond angles. Hydrogen does not attack the 
Bi nanoline, while ozone will oxidise around the nano- 
line at room temperature, leaving it intact. Likewise am- 
monia preferentially attacks the silicon, although some 
damage to the nanoline occurs at small fluxes. 



The determination of both the physical and elec- 
tronic structure of the nanoline has demonstrated the 
importance of a close interaction between experiment 
and theory. We note that all the recent publications on 
the Bi nanoline have been a collaboration between the 
two. The method of tightbinding has proved to be a valu- 
able tool for quickly searching through a large variety of 
possible structures, as in for instance the Haiku struc- 
ture, or for relaxing a large number of atoms, as in mod- 
elling of the nanoline kink energies. At the same time, 
the tightbinding method retains a quantum mechanical 
description of the energetics, so that the calculated en- 
ergies are in most cases in reasonable agreement with 
those calculated using DFT. The Bi nanoline has ben- 
efitted from an unusual degree of theoretical investiga- 
tion. A similar level of theoretical work for the rare-earth 
nanowires is likely to deepen greatly the understanding 
of the different nuclcation and growth phenomena that 
have been identified in that family of systems. 

With a significant body of work devoted to the fabri- 
cation and structure of the nanoline systems, attention 
should be paid to their possible applications. The length, 
straightness and perfection of self-assembled nanolines 
lend themselves to use as a nanoscale template. The most 
obvious candidate systems are the rare-earth nanowires, 
which are themselves metallic, so that an array of par- 
allel nanowires would readily form the contacts for a 
nanoelectronic device, perhaps using the "crossbar" ar- 
chitecture. The Er silicides appear also to act as pref- 
erential adsorption sites for metals such as Pt, so that 
the properties of these contacts is not limited to that of 
the nanowires themselves. However, the growth of these 
nanowires on a semiconducting substrate, with no indi- 
cation as yet that the substrate could subsequently be 
oxidised, might limit their application in nanoelectronics 
(for which a high-quality insulating substrate is crucial). 
For the Bi nanolines, oxidation of the surrounding silicon 
with ozone, followed by passivation of the surface with 
hydrogen, or nitridation using ammonia, is a plausible 
route for the fabrication of a single-nm ID template on 
an insulating surface, considerably smaller than any of 
the rare-earth nanowires, and well within the size range 
where quantum scale effects might be expected to occur. 
By the deposition of metals, active molecules, nanoparti- 
cles, or other species onto this template, the fabrication 
of nanowires, or arrays of nanoparticles with interest- 
ing optical, electronic or magnetic properties could be 
achieved, with atomic precision. This area is being ac- 
tively pursued by the authors. 

Acknowledgements We are happy to acknowledge useful dis- 
cussions with Chigusa Ohbuchi and Wataru Yashiro, as well 
as permission from various authors to reproduce figures as 
noted in the text. This study was performed through Special 
Coordination Funds for Promoting Science and Technology 
from the MEXT, Japan (ICYS and Active Atom- Wire Inter- 
connects). DRB is funded by the Royal Society. 



30 



J.H. G.Owen et al. 



References 

1. D.R.BOWLER J. Phys. Cond. Matt. 16 (2004) R721. 

2. J. V. Barth, G. Costantini and K. Kern Nature 437 
(2005) 671. 

3. W.Chen and H.Ahmed Appl. Phys. Lett. 62 (1987) 
1499. 

4. T. Hashizume, S. Heike, M. I. Lutwyche, 
S. Watanabe et al. Jpn. J. Appl. Phys. PL 2 35 
(1996) L1085. 

5. J. R. Heath, P.J.Kuekes, G.S. Snider and 
R.S.Williams Science 280 (1998) 1716. 

6. J.Nogami, S.L.Park and C.F.Quate Phys. Rev. B 36 
(1987) 6221. 

7. O. GuRLU, O. A. O. Adam, H. J. W. Zandvliet and 
B. Poelsema Appl. Phys. Lett. 83 (2003) 4610. 

8. M.Naitoii, H.Siiimaya, S.Nishigaki, N. Oishi et al. 
Surf. Set. 377-379 (1997) 899. 

9. K.MiKi, J. E.G. Owen, D.R. Bowler, G.A.D.Briggs 
et al. Surf. Sci. 421 (1999) 397. 

10. A. A. Baski, K. M. Saoud and K. M. Jones Appl. 
Surf Sci. 182 (2001) 216. 

11. C. Preinesberger, S. Vandre, T. Kalka and 
M. Dahne-Prietsch J. Phys. D: Appl. Phys. 31 (1998) 
L43. 

12. K.Sagisaka, D.Fujita and G.Kido Phys. Rev. Lett. 
91 (2003) 146103. 

13. O. KuBO, J. T. Ryu, H. Tani, T. Harada et al. 
Appl. Surf. Sci. 169 (1999) 93. 

14. J.T.Wang, H.Mizuseki, Y.Kawazoe, T. Hashizume 
et al. Phys. Rev. B 67 (2003) 193307. 

15. J. Wang, E.G.Wang, D.S.Wang, H.Mizueki et al. 
Phys. Rev. Lett. 94 (2005) 226103. 

16. C. J. Chen Introduction to Scanning Tunneling Mi- 
croscopy Oxford Scries in Optical and Imaging Sciences 
Oxford University Press (1993). 

17. R. Weisendanger Scanning Probe Microscopy and 
Spectroscopy: Methods and Applications Cambridge 
University Press (1994). 

18. A. R. Avery, H. T. Dobbs, D. M. Holmes, B. A. 
Joyce et al. Phys. Rev. Lett. 79 (1997) 3938. 

19. K.MiKi, Y.MoRiTA, H.Tokumoto, T.Sato et al. Ul- 
tramicroscopy 42-44 (1992) 851. 

20. B.S.Swartzentruber Phys. Rev. Lett. 76 (1996) 459. 

21. M. A. Rezaei, B. C. Stipe and W. Ho J. Chem. Phys. 
110 (1999) 3548. 

22. Y.Wang, M.J.Bronikowski and R.J.Hamers Surf. 
Sci. 311 (1994) 64. 

23. J.H. G.Owen, K.Miki, D.R. Bowler, C.M.Goringe 
ET AL. Surf Sci. 394 (1997) 79. 

24. B.Voigtlander, T.Weber, P.Smilauer and 
D.E.WOLF Phys. Rev. Lett. 78 (1997) 2164. 

25. J. H.G.Owen, K.Miki, D.R. Bowler, C.M.Goringe 
et al. Surf Sci. 394 (1997) 91. 

26. I.Goldfarb, p. Hayden, J. Owen and G. Briggs 
Phys. Rev. Lett. 78 (1997) 3959. 

27. R. M. Martin Electronic Structure: Basic Theory and 
Practical Methods Cambridge University Press (2004). 

28. M. W. FiNNis Interatomic Forces in condensed matter 
Oxford University Press (2003). 

29. A. SzABO and N. S. Ostlund Modern Quantum 
Chemistry McGraw-Hill (1989). 



30. C. M. Goringe, D. R. Bowler and E. Hernandez 
Rep. Prog. Phys. 60 (1997) 1447. 

31. P. Ordejon Comput. Mater. Sci. 12 (1998) 157. 

32. O. F. Sankey and D. J. Niklewski Phys. Rev. B 40 
(1989) 3979. 

33. A. P. Horsfield and A. M. Bratkovsky J. Phys.: 
Condens. Matter 12 (2000) Rl. 

34. A. P. Sutton, M. W. Finnis, D. G. Pettifor and 
Y. OllTA J. Phys. C 21 (1988) 35. 

35. D.R.BOWLER and J. E.G. Owen J. Phys. Cond. Matt. 
14 (2002) 6761. 

36. P. Hohenberg and W. Kohn Phys. Rev. 136 (1964) 
B862. 

37. W. Kohn and L. J. Sham Phys. Rev. 140 (1965) 
A1133. 

38. M. D. Segall, p. J. D. Lindan, M. J. Probert, C. J. 
Pickard et al. J. Phys. -.Condens. Matter 14 (2002) 
2717. 

39. S. Goedec'KER Rev. Mod. Phys. 71 (1999) 1085. 

40. X.-P. Li. R. W. Nunes and D. Vanderbilt Phys. Rev. 
B 47 (1993) 10891. 

41. S. Goedecker and L. Colombo Phys. Rev. Lett. 73 
(1994) 122. 

42. J. L. Fattebert and J. Bernholc Phys. Rev. B 62 
(2000) 1713. 

43. T. OZAKI and K. Terakura Phys. Rev. B 64 (2001) 
195126. 

44. J. M. SOLER, E. Artacho, J. D. Gale, A. Garcia 
et al. J. Phys.: Condens. Matter 14 (2002) 2745. 

45. D.R.BOWLER, T.Miyazaki and M.J.GlLLAN J. 
Phys.: Condens. Matter 14 (2002) 2781. 

46. C.-K. Skylaris, P. D. Haynes, A. A. Mostofi and 
M. C. Payne J. Chem. Phys. 122 (2005) 084119. 

47. G. A. D. Briggs and A. J. Fisher Surf Sci. Rep. 33 
(1999) 1. 

48. W. A. HOFER, A. S. Foster and A. L. Siiluger Rev. 
Mod. Phys. 75 (2003) 1287. 

49. J. Tersoff and D. R. Hamann Phys. Rev. B 31 (1985) 
805. 

50. B.S.Swartzentruber Surf. Sci. 386 (1997) 195. 

51. G. Brocks, P. J. Kelly and R. Car 70 (1993) 2786. 

52. M. M. R. Evans and J. Nogami Phys. Rev. B 59 
(1999) 7644. 

53. Z.-C. Dong, D. Fujita and H. Nejoh Phys. Rev. B 
63 (2001) 115402. 

54. N. Takeuchi Phys. Rev. B 63 (2001) 035311. 

55. G.P.Srivastava and R.H.Miwa Appl. Surf. Sci. 235 
(2004) 293. 

56. M. A. Albao, M. M. R. Evans, J. Nogami, D. Zorn 
ET AL. Phys. Rev. B 72 (2005) 035426. 

57. A. Kida, H. Kajiyama, S. Heike, T. Hashizume 
ET AL. Appl. Phys. Lett. 75 (1999) 540. 

58. J.-Z. Wang, J.-F. Jia, X. Liu, W.-D. Chen et al. 
Phys. Rev. B 65 (2002) 235303. 

59. J.-L. Li, X.-J. Liang, J.-F. Jia, X. Liu et al. Appl. 
Phys. Lett. 79 (2001) 2826. 

60. Y. P. Zhang, L. Yang, Y. H. Lai, G. Q. Xu et al. 
Appl. Phys. Lett. 84 (2004) 401. 

61. J.-L. Lin, D. Y. Petrovykh, J. Viernow, F. K. Men 
et al. J. Appl. Phys. 84 (1998) 255. 

62. R. Losio, K. N. Altmann, A. Kirakosian, J.-L. Lin 
ET AL. Phys. Rev. Lett. 86 (2001) 4632. 



Self-assembled nanowires on semiconductor surfaces 



31 



63. S. S. Lee, N. D. Kim, C. G. Hwang, H. J. Song 
ET AL. Phys. Rev. B 66 (2002) 115317. 

64. J. R. Ahn, H. Yeom, H. S. Yoon and I.-W. Lyo Phys. 
Rev. Lett. 91 (2003) 196403. 

65. J. N. Grain, A. Kirakosian, K. N. Altmann, 
G. Bromberger et al. Phys. Rev. Lett. 90 (2003) 
176805. 

66. A.Cho and J.M.Seo Surf. Sci. 565 (2004) 14. 

67. P. Segovia, D. Purdie, M. Hengsberger and 
Y. Baer Nature 402 (1999) 504. 

68. T.-C. Shen, C. Wang, G. G. Abeln, J. R. Tucker 
ET AL. Science 268 (1995) 1590. 

69. S. Watanabe, Y. a. Ono, T. Hasiiizume and 
Y. Wada Phys. Rev. B 54 (1996) R17308. 

70. S. Watanabe, Y. A. Ono, T. Hashizume and 
Y. Wada Surf. Sci. 386 (1997) 340. 

71. T. Hitosugi, T. Hasiiizume, S. Heike, S. Watanabe 
et al. Jpn. J. Appl. Phys. Pt. 2 36 (1997) L361. 

72. T. Hitosugi, S. Heike, T. Onogi, T. Hashizume 
ET AL. Phys. Rev. Lett. 82 (1999) 4034. 

73. C. F. Bird and D. R. Bowler Surf Sci. 531 (2003) 
L351. 

74. D. R. Bowler and A. J. Fisher Phys. Rev. B 63 
(2001) 035310. 

75. M. ToDOROVic, A. J. Fisher and D. R. Bowler J. 
Phys.-.Condens. Matter 14 (2002) L749. 

76. G. F. Bird, A. J. Fisher and D. R. Bowler Phys. 
Rev. B 68 (2003) 115318. 

77. D. P. Adams, T. M. Mayer and B. S. Swartzentru- 
BER Appl. Phys. Lett. 68 (1996) 2210. 

78. T. Hasiiizume, S. Heike, M. I. Lutwyche, 
S. Watanabe et al. Surf. Set. 386 (1997) 161. 

79. T.-C. Shen, G. Wang and J. R. Tucker Phys. Rev. 
Lett. 78 (1997) 1271. 

80. M. Sakurai, G. Thirstrup and M. AoNO Phys. Rev. 
B 62 (2000) 16167. 

81. M. G. Hersam, N. p. Guisinger and J. W. Lyding 
Nanotechnology 11 (2000) 70. 

82. G. C. Abeln, M. G. Hersam, D. S. Thompson, S.-T. 
Hwang et al. J. Vac. Sci,. Technol. B 16 (1998) 3874. 

83. G. G. Abeln, S. Y. Lee, J. W. Lyding, D. S. Thomp- 
son ET al. Appl. Phys. Lett. 70 (1997) 2747. 

84. G. P. LoPiNSKi, D. D. M. Wayner and R. A. 
WOLKOW Nature 406 (2000) 48. 

85. F. P. Netzer J. Phys.: Condens. Matter 7 (1995) 991. 

86. Y. Ghen, D. a. a. Ohlberg, G. Medeiros-Ribeiro, 
Y. A. Ghang et al. Appl. Phys. Lett. 76 (2000) 4004. 

87. Y. Ghen, D. A. A. Ohlbery, G. Medeiros-Ribeiro, 
Y. A. Ghang et al. App. Phys. A 75 (2002) 353. 

88. G. Ghen, J. Wan, J. Yang, X. Ding et al. Surf. Sci. 
513 (2002) 203. 

89. Y. Ghen, D. A. A. Ohlberg and R. S. Williams J. 
Appl. Phys. 91 (2002) 3213. 

90. R. Ragan, Y. Ghan, D. A. A. Ohlberg, 
G. Medeiros-Ribeiro et al. J. Cryst. Growth 
251 (2003) 657. 

91. L. Fitting, M. G. Zeman, W.-G. Yang and R. J. 
Nemanich J. Appl. Phys. 93 (2003) 4180. 

92. D. Lee, D. K. Lim, S. S. Bae, S. Kim et al. App. 
Phys. A 80 (2005) 1311. 

93. S. Harako, K. Kouno, S. Komuro, A. Ohata et al. 
J. Cryst. Growth 275 (2005) c2263. 



94. J. NoGAMi, B. Z. Liu, M. V. Katkov, G. Ohbuchi 
et al. Phys. Rev. B 63 (2001) 233305. 

95. G. Preinesberger, S. K. Becker, S. Vandre, 
T. Kalka et al. J, Appl. Phys. 91 (2002) 1695. 

96. B. Z. Liu and J. Nogami Nanotech. 14 (2003) 873. 

97. B. Z. Liu and J. Nogami Surf. Sci. 540 (2003) 136. 

98. Z. He, M. Stevens, D. .J. Smith and P. A. Bennett 
Appl. Phys. Lett. 83 (2003) 5292. 

99. B. Z. Liu and J. Nogami J. Appl. Phys. 93 (2003) 593. 

100. Z. He, D. J. Smith and P. A. Bennett Phys. Rev. 
Lett. 93 (2004) 256102. 

101. Z. He, D. J. Smith and P. A. Bennett Phys. Rev. B 
70 (2004) 241402(R). 

102. B. G. Harrison, P. Ryan and J. J. Boland Surf. Sci. 
582 (2005) 79. 

103. G. Ohbuchi and J. Nogami Phys. Rev. B 66 (2002) 
165323. 

104. G. Ohbuchi and J. Nogami Surf. Sci. 579 (2005) 157. 

105. M. Kuzmin, R. E. Perala, R.-L. Vaara, P. Laukka- 
nen et al. J. Cryst. Growth 262 (2004) 231. 

106. M. V. Katkov and J. Nogami Surf Sci. 524 (2003) 
129. 

107. M. Katkov and J. Nogami Bull. Am. Phys. Soc. 47 
(2002) 283. 

108. R. E. Perala, M. Kuzmin, P. Laukkanen, R.-L. 
Vaara et al. Surf Sci. 584 (2005) 8. 

109. M. Kuzmin, P. Laukkanen, R. E. Perala, R.-L. 
Vaara et al. Appl. Surf Sci. 222 (2004) 394. 

110. J. L. McGhesney, a. Kirakosian, R. Bennewitz, 
.J. N. Grain et al. Nanotechnology 13 (2002) 545. 

111. M. Tanaka, F. Ghu, M. Shimojo, M. Takeguchi 
et al. Appl. Phys. Lett. 86 (2005) 183104. 

112. J.-F. Lin, J. P. Bird, Z. He, P. A. Bennett et al. 
Appl. Phys. Lett. 85 (2004). 

113. M. Stevens, Z. He, D. J. Smith and P. A. Bennett 
J. Appl. Phys. 93 (2003) 5670. 

114. Z. He, M. Stevens, D. J. Smith and P. A. Bennett 
Surf Sci. 524 (2003) 148. 

115. H. Okino, L Matsuda, R. Hobara, Y. Hosomura 
et al. Appl. Phys. Lett. 86 (2005) 233108. 

116. G. L. Molnar, G. Peto, E. Zsoldos, N. Q. Khanh 
et al. Thin Solid Films 317 (1998) 417. 

117. B.G.Harrison and J. J. Boland Surf Sci. In proofs 
(2005). 

118. M. Kuzmin, R. E. Perala, P. Laukkanen, R.-L. 
Vaara et al. Appl. Surf Sci. 214 (2003) 196. 

119. B. Z. Liu and J. Nogami Surf Sci. 488 (2001) 399. 

120. R. Ragan, S. Kim, X. Li and R. S. Williams Appl. 
Phys. A 80 (2005) 1339. 

121. M. Kuzmin, R. E. Perala, P. Laukkanen and I. J. 
Vayrynen Phys. Rev. B 72 (2005) 085343. 

122. Z. He, D. J. Smith and P. A. Bennett Appl. Phys. 
Lett. 86 (2005) 143110. 

123. M. Fujimori, S. Heike, Y. Terada and 
T. Hashizume Nanotechnology 15 (2004) S333. 

124. D. Lee and S. Kim Appl. Phys. Lett. 82 (2003) 2619. 

125. K.MiKi, D.R. Bowler, J. H.G.Owen, G.A.D.Briggs 
et al. Phys. Rev. B 59 (1999) 14868. 

126. M.Naitoh, M. Takei, S. Nishigaki, N. Oishi et al. 
Jpn. J. Appl. Phys. Pt.l 39 (2000) 2793. 

127. J. H.G.Owen, K.Miki, H.Koh, H.W.Yeom et al. 
Phys. Rev. Lett. 88 (2002) 226104. 



32 



J.H.G.Owen et al. 



128. D.R. Bowler Phys. Rev. B 62 (2000) 7237. 

129. J.H.G.Owen, D.R.Bowler and K.Miki Surf. Sci. In 
press (2005). 

130. R. MiwA, T. Schmidt and G. Srivastava Surf. Sci. 
507-510 (2002) 368. 

131. R.H.MiwA and G.P.Srivastava Phys. Rev. B 66 
(2002) 235317. 

132. J.H.G.Owen, D.R.Bowler and K.Miki Surf Sci. 
Lett. 527 (2003) L177. 

133. J. MacLeod, R. Miwa, G. Srivastava and 
A. McLean Surf. Sci. 576 (2005) 116. 

134. M.Shimomura, K.Miki, T.Abukawa and S.Kono 
Surf Sci. Lett. 447 (2000) L169. 

135. J.H.G.Owen, D.R.Bowler and K.Miki Surf Sci. 
Lett. 499 (2002) L124. 

136. J. M. MacLeod and A. B. McLean Phys. Rev. B. 70 
(2004) R041306. 

137. R.H.MiwA, J. M. MacLeod, G.P.Srivastava and 
A. B. McLean Appl. Surf Sci. 244 (2005) 157. 

138. R. H. Miwa, J. M. MacLeod, A. B. McLean and 
G. P. Srivastava Nanotechnology 16 (2005). 

139. S.B.Zhang, W.E.McMaiion, J.M.Olson and S.-H. 
Wei Phys. Rev. Lett. 87 (2001) 166104. 

140. J.H.G.Owen, D.R.Bowler and K.Miki Mat. Sci. 
Technol. 20 (2004) 955. 

141. O. Sakata, W. Yashiro, D. R. Bowler, A. Kitano 
ET al. Phys. Rev. B 72 (2005) 121407(R). 

142. S.Yagi, W. Yashiro, K.Sakamoto and K.Miki Surf. 
Sci. In press (2005). 

143. H.Matsuiiata, K.Sakamoto and K.Miki J. El. Mi- 
crosc. 53 (2004) 325. 

144. A.Saito, K.Matora, T.Kurata, J.Maruyama 
ET al. Jpn. J. Appl. Phys. Pt.l 42 (2003) 2408. 

145. M.Naitoh, H.Shimaya, S.Nishigaki, N. Oishi et al. 
Appl. Surf. Sci. 142 (1999) 38. 

146. J. M. MacLeod, C. P. Lima, R. H. Miwa, G. P. 
Srivastava et al. Mat. Sci. Technol. 20 (2004) 951. 

147. H. J. W. Z. H. W. D. J. Wentink, A. van Silfhout 
and H. B. Elswijk Phys. Rev. Lett. 70 (1993) 2122. 

148. J.H.G.Owen, D.R.Bowler, C.M.Goringe, K.Miki 
ET al. Surf Sci. Lett. 341 (1995) L1042. 

149. M.Naitoh, M.Takei, S.Nishigaki, N. Oishi et al. 
Surf. Sci. 482-485 (2001) 1440. 

150. S. R. SciiOFiELD, N. J. Curson, J. L. O'Brien, M. Y. 
Simmons et al. Phys. Rev. B 69 (2004) 085312. 

151. A. D. Becke and K. E. Edgecombe J. Chem. Phys. 
92 (1990) 5397. 

152. D. R. Bowler Surf. Sci. (2005) in press. 

153. J.H.G.Owen and K.Miki (2005) In Preparation. 

154. K.Nakamura, S.Ichimura, A.Kurokawa, K.Koike 
ET al. J. Vac. Sci. Technol. A 17 (1999) 1275. 

155. G.-M. Rignanese and A. Pasquarello Appl. Phys. 
Lett. 76 (2000) 553. 

156. I.F.KovAL, P.V.Melnik, N.G.Nakiiodkin, M. Py- 
atinsky et al. Surf. Sci. Lett. 384 (1997) L844. 

157. T.V.Afanasieva, I.F.Koval and N.G.Nakhodkin 
Surf Sci. 507-510 (2002) 789. 

158. J.H.G.Owen, D.R.Bowler, S.Kusano and K.Miki 
Phys. Rev. B 72 (2005) 113304. 

159. T. Itoh, S. Kashirajima, M. Naitoh, S. Nishigaki 
et al. Appl. Surf Sci. 244 (2005) 161. 

160. A. V. ZoTOv, A. A. Saranin, K. V. Ignatovicii, 
V. G. LiFSHlTS ET AL. Surf. Sci. 391 (1997) L1188. 



161. K. MiKi, H.Matsuiiata, K. Sakamoto, G. A. D. 
Briggs et al. Microscopy of Semiconducting Materials 
XL Lnst. Phys. Conf Ser. 164 (1999) 167. 

162. O. Sakata, M. Takata, H. Suematsu, A. Matsuda 
et al. Appl. Phys. Lett. 84 (2004) 4239. 



